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Abstract. This article presents the first computation of the complete 
bispectrum of the cosmic microwave background temperature anisotropics arising 
from the evolution of all cosmic fluids up to second order, including neutrinos. 
Gravitational couplings, electron density fluctuations and the second order 
Boltzmann equation are fully taken into account. Comparison to limiting cases 
that appeared previously in the literature are provided. These are regimes for 
which analytical insights can be given. The final results are expressed in terms 
of equivalent fj^^ for different configurations. It is found that for moments up to 
^max = 2000, the signal generated by non-linear effects is equivalent to fj^^ ~ 5 
for both local-type and cquilateral-type primordial non-Gaussianity. 



1. Introduction 

The Cosmic Microwave Background (CMB) anisotropics arc now observed with a 
high precision and have become a key observation of modern cosmology. They are in 
particular very precious to constrain the theories of the primordial Universe [1]. So 
far the temperature anisotropics have been found to have statistical properties that 
are compatible with Gaussian statistics [1]. The CMB data can therefore entirely be 
captured in its power spectrum and the latter has been used to set constraints on 
the cosmological parameters and on the shape of the inflationary potential. There is 
hope however that future observations such as Planck [2] , that will provide better data 
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on the statistical properties of the temperature and polarization fields, open a new 
window on the physics of the early Universe with the use of higher-order statistical 
properties of the CMB sky such as its bispectrum. 

In the analysis of those data, one should keep in mind though that the properties 
of the CMB temperature and polarization anisotropies depend both on the properties 
of the initial conditions and on their evolution. As long as measurements are restricted 
to second-order statistics such as the angular correlation function or the power spectra, 
a linear perturbation theory suffices, for the required precision, to relate the (2- 
dimensional) angular power spectrum to the (3-dimensional) initial power spectrum 
of the metric perturbations at the end of inflation. It has thus been understood 
that the characteristic features observed in the CMB temperature spectrum originate 
from the developments of acoustic oscillations encoded in the linear transfer [3J H] 
function while its overall amplitude and its scale dependence are fixed by the initial 
power spectrum, the shape of which agrees with the predictions from an inflationary 
era [SJ [B] . At this level of description, all aspects are now fully understood and is 
part of textbooks [71 HI [9l [TOl [Til HI] and since the metric and matter perturbations 
are linearized, any model that predicts Gaussian initial conditions, as standard single 
field inflation, is expected to produce Gaussian statistical properties for the CMB 
temperature field. 

In general bispectra arise whenever non-linear mode couplings are at play during 
the cosmological evolution and general relativity being in essence a non-linear field 
theory, deviations from Gaussianity are expected to be ubiquitous, arising either 
from the inflationary era (and thus called primordial no n- Gaussianity) or from the 
post-inflationary evolution. Gravity mediated couplings are however generally small 
and in full agreement with the current data that clearly favors only mild non- 
Gaussianities if any [1]. In particular, it is now widely accepted that standard 
single field inflation cannot produce signiflcant non-Gaussianities since its amplitude 
is mostly dictated by gravity induced couplings [13j . On the other hand, significant 
deviations from Gaussianity can arise from non-standard kinetic terms, for which 
models based on the Dirac-Born-Infeld action are typical models [Tl] or in the 
context of of multiple-field inflation specially when non-gravity type couplings are 
at play [15l[T6l[T7l[T8l[19l[20l[2ll[22|. It is therefore generally admitted that non- 
Gaussianity searches can open a window on the details of the inflationary mechanism 
at work in the early Universe (see Refs. [23l[2ll[25] for general reviews). 

The level with which primordial non-Gaussianity could actually be detected is 
however still largely debated, mainly because this source of non-Gaussianity is in 
"competition" with the couplings induced by non-linear effects throughout the whole 
recombination and photon propagation processes. This has motivated a series of 
general studies aiming at characterizing the bispectrum [26l [27l EH [29l [30l [HI [32j 
[33l[Hl[35l[36l[3ll[38l[3i[40l[4ll[42] and even the trispectrum [HI [43l [44] to be 
expected from the observation made by CMB experiments. While the identification of 
the mode couplings (at the quantum level) during the inflationary phase has been set 
on secure grounds [131 [IB [Ml [13 [48] , the evolution of the cosmological perturbations 
during the post-inflationary era is still largely unexplored, the primordial statistical 
properties being often related to the observed statistical properties only through a 
linear transfer function. However, in order to relate the angular bispectrum of the 
observed CMB to the spatial bispectrum of the metric perturbations at the end of 
inflation, one needs to derive transfer functions up to second order that incorporate 
all type of couplings. 
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Effect 


Linear evolution 


Non-linear evolution 


Primordial 


Primordial power spectrmn P{k) 


Primordial /^j^ 


Primary (early) 


Sachs- Wolfe and Doppler effects 


This article 


Secondary (late) 


Integrated Sachs- Wolfe effect, 
reionization 


Lcnsing-ISW correlation 



Table 1. Classification of the linear and non-linear effects. 



The work we present in this article is the end result of a task which was initiated 
in Ref. [49] , that was followed by some partial reanalysis [35l \EU[ [51] , In this article 
we described what we thought is the main mechanism at play for the generation 
of the bispectra at small angular scales and for which we could give some physical 
insights. The present article extends this analysis using the theoretical developments 
of the non- linear perturbation theory described in Refs. [52l [53l [54] , and in particular 
concerning the second order Boltzmann equation that we need in order to describe 
the evolution of radiation and neutrinos. Its goal is three-fold: (i) to support the 
predictions and analytic understanding obtained on small scales in Ref. [49] by a 
full numerical integration of the second order Boltzmann equation, without neglecting 
any term, (ii) to compute the expected bispectrum on intermediate scales, and (iii) 
to revisit the large-scale behaviour of the bispectrum which has been presented in 
the previous literature [32l [40] [55]. We emphasize that the amplitude of this non- 
Gaussianity is completely fixed once the amplitude and power spectrum of the first 
order scalar perturbations are constrained so that this offers a definite prediction of 
the minimal amount of non-Gaussianity expected in any CMB observation, that is 
provided Einstein theory remains a good description of gravity. 

Let us mention that our present analysis refines our previous descriptions by 
classifying the evolution effects, linear or not, in two categories. First, the primary or 
early effects arise from the evolution of the perturbations in the radiation dominated 
era after horizon crossing, during recombination between protons and electrons and 
until the potentials are constant in their linear evolution. Then secondary or late 
effects arise at late time when the potential starts to evolve due the late time 
acceleration of the Universe or due to reionisation effects (see Ref. [56] and references 
therein). This separation is unambiguous for the linear evolution and is therefore 
pertaining to the power spectrum computation but also to the calculation of the 
bispectrum that mix linear and second order source terms. This is actually clear 
from Fig. [1] below, and it is summarized in table [TJ This present work focuses 
on the non-linear evolution of the field during the early period and the resulting 
bispectra it induces. In particular we compute the shape, amplitude and bispectrum 
of the temperature anisotropietl^ due to these effects and eventually compare them to 
secondary sources, and specifically the ISW-lensing effects, and to primordial coupling 
effects (through standard parameterizations). In this work the non-linear evolution of 
the fluids, including neutrinos, is therefore treated exactly up to second-order in the 
linear perturbations until reionization is complete. That includes of course the use of 
the second order Boltzmann equation. Adiabatic initial conditions are also assumed. 
All the results described here can actually be reproduced from a Mathematica code 

I Note that these results actually depend on the actual definition of the temperature one uses (see 
the discussion of S 14.11 1. This is due to the spectral distortion that second order effects necessarily 
induce as described in Ref. 1571 . 
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which is freely available - with its documentation - on the webpage |58] . 

The article is organized as follows. Section 2 describes the main concepts of 
the non-linear cosmological perturbation theory while all equations arc gathered 
in [Appendix A Then, in Section 3 we determine the initial conditions for both the 
metric variables and cosmic fluids (baryon, cold dark matter, photons, neutrinos). 
In Section 4 we describe, after a careful definition of the temperature in § 14.11 the 
numerical integration based on the flat sky approximation presented in Section 3. 
The numerical results concerning the bispectrum are discussed in Section 5 while in 
Section 6 we provide an analytical understanding of these results in various limiting 
cases. 



2. Non-linear Perturbation Theory 

We follow Ref. [53] for the notation. This section summarizes our main choices and 
conventions and we refer the reader to the latter reference for any further details. We 
first detail the description of the metric in § 12.11 of the matter fields in § 12.21 and then 
describe the structure of the perturbation equations in § 12.31 



2.1. Metric perturbations 

At lowest order, we assume that the Universe is well described by a Friedmann- 
Lemaitre (FL) space-time with Euclidean spatial sections and with scale factor a. The 
Universe is then described by a perturbed space-time around this FL background. In 
the Newtonian gauge (often also named Poisson gauge), the form of the perturbed 
metric is then given by 

ds^ = .g^^d-T^dx"^ (2.1) 
= a{r]f{ -{1 + 2$)d772 + 2Bidx'di] + [(1 - 2^)Sij + 2Hij]dx' dx''} , 

with d^Bj ~ d^Hjj = H^j = 0. Here, indices of the type I,J,K ... are spatial 
coordinates indices running from 1 to 3, and rj is the conformal time. We also define 
the comoving Hubble factor % = a' /a where a prime denotes a derivative with respect 
to 77. Each of the perturbation variables has an order by order expansion of the 
type [Sn] 

X = X + X'^^^ + ]^X^^K (2.2) 

In the case of metric perturbations, the definition is made such that their background 
value vanishes. Additionally, since vector modes fail to be generated at first order in 
the standard models of inflation [3 O [10] , we only consider the second order vector 
modes Bi^'^\ Similarly, the amplitude of first order tensor modes (i.e. gravitational 

waves) is expected to be sufficiently low [50] so that we can safely neglect j , and 

(2) 

we will consider only Hj j . 



2.2. Describing the matter content 

In the concordance model of cosmology, the matter content of the Universe includes 
relativistic particles (radiation (photons), neutrinos) and non relativistic particulcs, 
which fall either in the category of baryons or cold dark matter. Once the equation 
of state for each species is specified, the background energy density is sufficient to 
characterize its state. They are denoted p^, pr, Pc and pb for neutrinos, radiation. 
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cold dark matter and baryons, respectively. Usually the (perfect) fluids contained in 
the Universe are characterized in the perturbed Universe by their energy density but 
also by their velocity. However this is not possible for relativistic particules and, in 
order to describe photons and neutrinos, we need to rely on a kinetic description based 
on the moments of their distribution function. We will thus describe the fluids as a 
special case of this statistical description when only the two first moments, which are 
related to the energy density and the velocity [53], are needed. 

2.2.1. The moments of species In general the moments of the distribution function 
are defined according to the spatial part of a local Minkowski space-time which itself is 
defined at any point of the space-time manifold by the value of a tetrad field. They are 
also taken with respect to an azimuthal direction that we arbitrarily align with a given 
Fourier mode when working out their spatial dependence in Fourier space. Wc thus 
describe the radiation with the brightness moments i^^^\\^,r|) and I™'"^'' (k, r/), but 
also with the moments describing its polarization which are separated into the electric- 
type moments E'^^^^ (k, 77) and f (k, rf), and the magnetic type moments B^^"^^ (k, 77) 
(which appear only at second order since we have assumed that the first order vector 
and tensor modes can be neglected). The polarization of radiation is generated through 
the interactions with baryons while the neutrinos, being collisionless, develop no 
polarization. As a consequence, we only need the brightness moments A/"™^^^ (k, 77) 

and J\f^^^^ (k, r/) to describe them (we also assume that the three species of neutrinos 
are massless and we describe them altogether though their mass splitting can have 
some interesting cosmological implications |611 162] ). 

The matter species, i.e. baryons and cold dark matter, are described by 
perfect fluids. The perturbation of their energy densities is defined according to the 
decomposition (|2.2|) . The spatial moments of the velocities of these fluids are noted 
v,n with TO = —1,0,1 and their perturbation is also defined according to Eq. (|2.2p with 
no background value thanks to the symmetries of the background space-time, which 
imply that all fluids have the four-velocity = since no preferred direction exists. 
Furthermore, since the first order vector modes are not generated at first order, we 
have v'^^l = 0. 

Additionaly, it will prove more convenient to work also with the moments of the 

(2) (2) 

vector and tensor perturbations, that is with and H)_{ instead of Bi and Hij. 
Again, we refer to Rcf. |53] for the exact definitions which are inspired from Rcf. |63j . 

2.2.2. Stress- energy tensor For a perfect fiuid, the stress-energy tensor is given by 

= pUf^u^ + P (.9^1, + u^Ui,) , (2.3) 

where the pressure is related to the energy density by an equation of state P = wp. We 
can see on this expression of the stress-energy tensor that we neglect the anisotropic 
stress and it thus applies for baryons and cold dark matter only. For radiation and 
neutrinos, the stress-energy tensor can be constructed from the distribution function 
and it is thus possible to relate its expression to the moments of the brightness (see 
Ref. [33] for instance). 

2.3. Field equations 

2.3.1. Structure of the equations Three equations are essential to understand the 
evolution of perturbations of the radiation. First, we need the Einstein equation 
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G^i, — nT^jj with K = SttG, in order to relate the perturbations of the metric to 
the ones of the different matter species. Then, the evolution of the distribution 
functions is dictated by the Boltzmann equation from which one can extract the 
evolution hierarchy of its moments. Finally, since the collision rate between photons 
and baryons is proportional to the fraction of free electrons, Xo, which is governing the 
recombination process, we need to derive the evolution equation of Xc up to the first 
order in the perturbations. These three sets of equations are detailed in [Appendix A| 
At linear order, the evolution equations are partial differential equations which 
are by construction linear in the perturbation variables. They formally take the form 

V,[xf\=Q, (2.4) 

where the X^^^ is the set of first order perturbation variables and Vi the set of 
differential operators acting on them. At second order the purely second order 
perturbation variables satisfy the same differential equations, but supplemented with 
source terms which are quadratic in the first order variables, i.e. they are of the form 

V,[xf^]=S,[xf\ (2.5) 

where the terms Si are operators (possibly differential) quadratic in their arguments. It 
is thus sufficient to write down the second order equations, that is the linear operators 
TDj and the source terms 5,;, in order to fully specify both the first and the second 
order equations. 

2.3.2. Structure of the solutions At first order, the space of solutions is spanned by 
the linear combinations of independent solutions since the operators "Di are linear. 
Then, the choice of the initial conditions selects a particular solution in this vector 
space of solutions. However, at second order, the space of solutions is spanned by the 
linear combinations of the same independent solutions to which a particular solution 
has to be added, thanks to the non vanishing source terms Si which break the linearity, 
so that it belongs to an affine space instead of a vector space. The choice of this 
particular solution is arbitrary so that we can decide to take the one which vanishes 
deep in the radiation dominated era. This particular solution will then be physically 
interpreted as the cumulative effect of the sources and this is actually what we mean 
when we compare the effects of non-linear evolution to those of the linear evolution, 
the latter being the solutions of the homogeneous equation. As explained in the 
introduction, the goal of this article is to assess the importance of the non-linear 
evolution effects on the temperature bispectrum with respect to the linear evolution 
effects which are now well under control from a theoretical and numerical point of 
view. 

2.3.3. Fourier space and convention From a technical point of view, we will solve 
numerically the evolution equations in Fourier space. We define the Fourier transform 
of any function X(x, ry) as 

^(k,^) = j ^^^X(x,,y)e-"'- (2.6) 

and we will use the expression of the Fourier transform of a product as the convolution 
[Xy](k) =/C{X(ki)r(k2)} , (2.7) 
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hence introducing the short hand notation 
where is the 3-diniensional Dirac distribution. 



(2.8) 



3. Initial conditions 



Since both radiation and neutrinos arc rclativistic, the ratio of their energy densities 
is conserved during the cosmological evolution, and we define the fraction of radiation 
among the rclativistic species by = Pr/iPr + Pu)- Similarly, cold dark matter and 
baryons being pressureless and coupled to gravity universally [53], the ratio of their 
energy densities is also conserved. We also define 



4 p. 



(3.1) 



In this section focusing on the initial conditions, all expressions are meant to hold only 
at an initial time r?init- We thus omit to specify it in order to simplify the notation. 
We now give the initial conditions for the matter perturbation variable in § 13.11 and 
the geometry perturbation variables in § 13. 2| both at first and second orders. 



3.1. Initial conditions for the matter variables 

We only consider in this article adiabatic initial conditions since isocurvature modes 
are highly constrained observationnally j65] but also because the single field models 
of inflation [51 [TU] predict adiabatic initial conditions. The link between isocurvature 
modes and multiple- field models remains to be investigated in depth [TU [66] . 

For adiabatic initial conditions, the velocity of all matter species is negligible both 
at the first and at second orders. It is indeed suppressed by a factor krji^it compared 
to 3>init- We thus take the initial condition 



,(1) 



,(2) 



0, with 



-1,0,1 



for baryons and cold dark matter and we take 



^(2) 



'(2) 



0, 



for radiation and neutrinos. Adiabaticity also implies that the ratios between 



(3.2) 



(3.3) 

1/3 



pV^ ^ pV'^ and remain constant. We deduce the following initial conditionf[§| 



3 pc 

Ip^ 
3 Pc 



3 Pb 
3 Pb 



rO(l) 



0(2) 



rO(l) 



(3.4a) 



(3.4&) 



Tr\k)=Nr\k). 



(3.4c) 
(3.4rf) 



§ Xq is not exactly the energy density since the frame used to define the multipoles is not necessarily 
comoving with the species. However the difference are quadratic in the velocity field, mainly because 
it arises from boost factors, and it is thus not important for the initial conditions. See Refs. 1531 1671 
for the link between these two quantities. 
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Now, using the Poisson equation (|A.ip deep in the radiation era, we can relate these 
initial conditions to those of the metric perturbations as 

-7-0(1) _ _ 2$(i) 



■-0 

7-0(2) 



1)2 



(3.5a) 
(3.56) 

2U1I does not vanish and this 



Additionally, the presence of neutrinos implies that ^ 
term is required in order to determine the initial conditions for the perturbations of 
the two gravitational potentials since $ — ^I* is proportional to A/'^ see Eq. (jA.3p . 
[Appendix C| details the initial conditions satisfied by the quadrupoles I^''^^ -^2 ^^^^ 



A/'^ while we always have 21 



■0(1) 

2 



3.2. Initial conditions for the metric variables 

As previously, we fix the initial conditions deep in the radiation dominated era for 
super- Hubble modes, so that we can expand all our equations in terms of fcr/init ^ 1. 
For adiabatic perturbations, we recall that on super-Hubble scales the comoving 
curvature perturbation is conserved. Its expressions at first and second orders 
are [Si [Si 

2 



7^(1) =^(1) 

7^(2) ^ $(2) 



2 

3{l + w)H 



3{l + w)p 



^'(i)+7^$(i) 



V]/'(2) ^ -^$(2) _ 4-^(j)(l)2 

^ 4 p^)^,. 



^(l)'2 

H 



(3.6a) 
(3.66) 



3(1 + w) p 

where w and are the equation of state and the adiabatic speed of sound of the 
total fluid. Deep in the radiation era the properties of the total fluid are close to the 
radiation+neutrinos fluid since pm ^ Pr ^-nd w = = 1/3. At first order, we obtain 
the initial conditions [7D] 



1 + ^(1-,/,-) 



$(1) = ^7^(l 



vl/(l)_cl>(l) = f(l_/jcl>(l), 

5 

which in particular implies 72. = ^'^^^ + <I>(i)/2. At second order, we obtain 

-5* + 5«'"' 



$(2) 


2 




H2$(i 




$(2) 


-I 


^^$(1)2^ 


h 2<I>(^) 















17^(2) 



where 



S3 - -n^A- 



fX,^'^ + {l-l.)N^ 



(3.7a) 
(3.76) 

(3.8a) 
(3.86) 

(3.9) 



[[ Note that the quadrupole moments, that is the moments with £ = 2 (J\f^ and X™) are different 
from the anisotropic stress of neutrinos and photons 1671 . In fact they match only if the frame in 
which the moments are taken is ahgned with the velocity of the species considered. 
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It can also be recast in the following set of two initial conditions 



2 f 3 



V[/(2) _ (j) 



(2) 



(3.10) 



(3.11a) 
(3.11&) 



In order to make link with previous works |26| . note that in the case where there are 
no neutrinos (/, = 1) we obtain 



$(1) ^ vi/d) = ±TZ^^) , 
3 

^(2) ^ ^^(2) ^ ^(1)2 _ ^-1 \^J^^^M^'¥^'> 

3 

$(2) = ^(2) + 4,j(i)2 _^ 3^-1 Idjx^jWgi^jW 



(3.12a) 
(3.126) 
(3.12c) 



For the vector and tensor modes, we have argued that they are respectively purely 
decaying modes or of negligible amplitude at first order, so that 

$S - i/iV = 0. (3.13) 

At second order, vector and tensor type perturbations exclusively arise from the 
integrated contribution of the quadratic sources so that they also vanish at sufficiently 
early times. We thus have 

= H^±l = 0. (3.14) 



3.3. Primordial non-Gaussianity and f^-^ 

In the previous section, we have not specified the value of TlP'\ Indeed this is what 
we would like to constrain by measuring the level of non-Gaussianity since different 
models of inflation have different predictions. In general, theoretical predictions are 
given in terms of the spectrum and of the bispectrum of TZ related to the 2- and 3-point 
functions by 

(7^(kl)7^(k2)) = (7^(l^(kl)7^(l)(k2)) = 5l (ki + ka) Pn{ki) , (3.15) 

and 

(7^(kl)7^(k2)7^(k3)) = i(7^(l)(kl)7^(l)(k2)7^(2)(k3)) + sym. 

= 51 (ki + k2 + k3) ff^ ^^(ki, k2, k3) , (3.16) 

where P-R.{k) is the primordial power spectrum and i^(ki,k2,k3) is a function which 
specifies the shape of the primordial non-Gaussianity and is its amplitude. 
Each possible type of primordial non-Gaussianity is characterized by a function 
i^(ki,k2,k3), and these functions are normalized such that they all have the same 
value in the Fourier configuration ki ^ k2 ~ k^. 
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Different limiting cases can be found in the literature |71) . the main one being 
the non-Gaussianity of the local type for which 

^^ioc(ki,k2,k3) = 2P7e(fci)PK(fc2) + 2sym. (3.17) 

In this case, it is equivalent to specify that 

7e(2)(k) = 2/c {/^^7^(l)(kl)7^(l)(k2)} (3.18) 

in Fourier space. In particular, it has been shown [13] that for slow-rolling single field 
inflationary models 

7^(2) ^ _27^(l)^ (3.19) 

that is ~ —1. Note that this is equivalent to impose that ~ with 

exp(— 27?,m) = 1 — 271 [13]. This is the initial condition that we shall assume in this 
work since we want to study the effects of evolution in the pursuit of non-Gaussianity 
and want to consider the minimum amount of primordial non-Gaussianity (see however 
the discussion of S 13.41 below). 

Note also that most expressions are given in terms of the power-spectrum of the 
gravitational potential 3>'^^' rather than 7?,*^^-'. The two power spectra are easily related 
thanks to 

Similarly, the primordial non-Gaussianity is more often expressed in terms of a /^^j^ 
parameter which is related to $ rather than TZ. In order to be consistent with previous 
literature in which neutrinos were not included at second order, we define f * through 

$(2)(k) = -2/C |/* ^ g (^1 + ± (1 - /J^] $(i)(ki)<i>(i'(k2)} , (3.21) 

which is inspired from Eq. p. flap when quadratic terms are ignored, that is when 
the relation between and 7?.'^-' is approximated to be linear, which is a valid 
approximation for large primordial non-Gaussinities. This is the non-Gaussianity 
parameter that we will discuss in the rest of this article. Historically, the definition 
p.2ip was made for = 1, and since the matter dominated era potential is related 
to the primordial one by $|j[at = 'Hj'^init large scales, then we would obtain 



mat 



(k) = -2/c {/* <I>W,(ki)$L'it(k2)} . (3.22) 



It is useful to notice that when we consider the neutrinos (/^, 7^ 1) for the concordance 
model [U 

3 

5 



~ 0.99 f$(i)V , (3.23) 



and thus, the same type of relation holds when considering neutrinos, but deep in the 
radiation dominated era. 

Another interesting family of non-Gaussian initial conditions, motivated by 
models such DBI inflation, is of the equilateral type where the function F peaks when 
ki = k2 = ks- The expression for such type of non-Gaussianity is given in Ref. [7^ 
and reads 

Fequi(ki,k2,k3) ^ 6 [-Pn{ki)Pn{k2) + 2 sym. - 2[Pn{ki)Pn{k2)Pnik3)f^ 

+P^\ki)P^\k2)Pn(k3) + 5sym.} . (3.24) 
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Similarly to what has been done in Eq. p.21[) we can also define a for the equilateral 
configuration. The local and equilateral type non-Gaussianities are interesting from a 
technical point a view, since there analytical expression in function of ki, ^2 and 
is factorizable in powers of these quantities. This implies that the optimal estimator, 
that we need to build to constrain this type of primordial non-Gaussianity, is faster 
than for a general shape of the primordial non-Gaussianity |31] . Even if the true 
primordial non-Gaussianity is not strictly of the local or equilateral type, it can be 
approximated by these ideal cases in order to produce constraints on the primordial 
non-Gaussianity in a fast and efficient manner j42j . 

3.4- Why we can concentrate on = — 1 and still be general 

Let us now stress that as long as one considers the bispectrum, the result is linear 
in the second order solution. This linearity shall not be confused with the non- 
linearity of the perturbation equations. Indeed, as long as the initial conditions at 
first order are identical, the second order initial conditions (j3.4bl I3.5M I3.12bl I3.12cp 
for the matter and metric perturbations are linear in TZ^^\ Since the solutions of the 
second order perturbation equation belong to an affine space, we can always write 
fT^ = fT^ + -.^ , ,. where the first quantity refers to the extra primordial non- 

•'nl ■'NLprim ■'NLlficld' M. J f 

Gaussianity compared to the standard one- field inflation prediction /^j^gpy It follows 
that the initial condition for TZ^^^ satisfies 

7^(^)(k) = 2/c {4^,p^.^„^7^(l)(kl)7^(l)(k2)} + 2/c {C^^^^^'H^^WH^^)} , 

that is TZ^'^^ (k) = TZ''pJ^^(k)+Tz'^l evolution C')- purely second order variables will 

thus be decomposed as X^'^'>{h) — x'^2^-^{'k) + X^^^ evolution C^) ^'^'^ ^^"^ final bispectrum 
will simply be the superposition of a primordial bispectrum and of the bispectrum 
induced by the non-linear evolution that is computed in this article. 

4. Computation of the spectrum and bispectrum 

We adopt the line of sight approach |63l [70l [73l [111 [75] in the resolution of the 
Boltzmann equation. This consists in considering that the observed brightness 
in a given direction of the CMB is the sum of the brightness of the emitting 
sources. It simplifies the numerical implementation since, at least at the first order in 
perturbations, only lower multipoles contribute to the sources. 

In this section, we first discuss the definition of the temperaturre in § 14.11 and 
then describe the expression of the emitting sources in § 14.21 (after having defined our 
choice of definition for the temperature in § 14. ip . Then we use a flat sky approximation 
(to be distinguished from the Limber approximation which is an extra approximation, 
which can be made or not) in order to simplify the numerical integration both for the 
spectrum in § l4.2l and for the bispectrum in § 14.51 We refer to Ref. [77] for the detailed 
derivation of the expressions summarized below and for a discussion of their accuracy. 

Defining the temperature 

The primary physical quantity related to the distribution function is the brightness, 
and we need to define the CMB temperature from this quantity. First we define the 
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brightness temperature fractional perturbation by 



4 \I 

where I = pj is the background value of the brightness, from which we can define the 
brightness temperature by Tx = T(l + G). This definition is inspired from the more 
general bolometric temperature T which is defined by 



, - (4.2) 

TJ I ^ ' 

The brightness temperature and the bolometric temperature agree at the background 
and first order level. Since the Boltzmann equation is derived in function of the 
brightness, the brightness temperature is the most straightforward quantity to use, 
but the bolometric temperature carries more meaning since it is the temperature of 
the black-body distribution which would have the same energy density as the actual 
distribution. However, it has been shown |57| that the non-linear dynamics sources 
a y-type spectral distortion and this would affect the brightness and thus both the 
brightness and the bolometric temperatures. Following our previous analysis j57j . 
it would be more natural to use the occupation number temperature defined as the 
temperature of a black-body which would have the same number density of photons 
as the observed distribution (to be contrasted with the bolometric temperature which 
is the temperature of the black body which carries the same energy density as the 
observed distribution). It is thus defined by 

Tn ^ —^—T ■ (4.3) 

(l + 42/)4 

Defined like this, this temperature is not affected by spectral distortions since Compton 
scattering conserves the number of photons in collisions, and it has a non-ambiguous 
signification with the type of distortions induced by the non-linear collision term |57j . 

These definitions agree at the background level and at first order in the 
perturbation so that Cj^"^^ = Cf^ = Cj"^" = T^Cf ®. Since they differ at second 
order in perturbations, we can define three different bispectra, depending on the choice 
of definition for the temperature. As we shall see in § 15. 4[ since at lowest order in the 
spectral distortions [57] 

T = TN{l + y), (4.4) 

the three bispectra are related to one another so that one only needs to compute one 
of them. In order for our results to be easily compared with the existing literature, we 
will use the bolometric temperature, although the occupation number temperature is, 
at least from a theoretical point of view, more natural to use. 

4-2. Emitting sources and transfer functions 

Given the definition of the multipoles used in Ref. [53] . the brightness temperature 
emitted by a multipole I™ (k, 77) is given by 

^QTmriKv), (4.5) 

where we have introduced 

QT{k)^Ni'Yr{k), (4.6) 
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k being the unit vector in the direction of k and where 



N.^i'f-^. (47) 

It follows that the total expression of the emitting sources is given by 

Wik,i^) EE exp(-f)^5r(k,^)Qr(k), (4.8) 

where f is the optical thickness which satisfies the integral equation 



Vo 



dr/f'(ry'), (4.9) 



that is df/d?7 = — f', and where the emitting sources multipoles SJ^ are defined 
in [Appendix B[ This expression is valid both at first and at second order, depending 
on whether we consider the first order or the second order sources. Since we want to 
determine how the primordial perturbations of the gravitational potential $ translate 
into the observed temperature, we are rather interested in the transfer functions 
defined by 

u;(i)(k,77)$(k,77i„it) = VF(i)(k,77), (4.10a) 

/C {u;(2)(ki, k2, r,)$(i) (ki, 77init)$^'' (ka, ?7init)} = W^^'> (k, r,) . (4.10&) 

As we shall see, w^^^(k, 77) and U''^^^(ki,k2,77) will be computed numerically. In 
principle the I™(k) are complex functions, but numerically we prefer to manipulate 
real numbers and real differential equations. At first order this is not an issue since only 
the scalar perturbations are non vanishing and 2^''^'' (k) remains unchanged under a 
rotation of the coordinates system around k, that is about the azimuthal direction used 
for the spherical harmonics. However, at second order, this is no more the case since 
for a given £ all the allowed m are to be taken into account. In order to perform the 
integration with real numbers, we thus choose the orientations of ki and k2 such that 
y^ki = fy^ko = 0; that is kiy = = 0. When the azimuthal direction of the spherical 
harmonics is k, ^^'"'^''(ki) and Z™^^''(k2) are then computed according to Eq. (jA.37|) . 
which ensures that they remain real. In the computation of the bispectrum, we will 
then need a general orientation of ki and k2. It can be obtained by rotating this 
configuration around the azimuthal direction k by an angle a(ki, k2). Defining I™^^^ 
in a similar way as Eq. (j4.10fep . 

/C [if^^ (ki , k2, ry)<i>(i) (ki , ?7i„it)$(^) (k2, ?7i„it) } = if^^ (k, r?) , (4.11) 
it can be related to the result of the (real) numerical integration Z^"^^j(ki, k2) by 

ir(ki,k2) = e'""(''-''=)ll^.,,i(ki,k2) . (4.12) 
Our notation a refers to the Euler angles. 
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4.3. Flat sky approximation 

In the flat sky approximation, the sky is expanded around a reference hne of sight 
with direction f . Any Fourier mode k can then be spht into a component fcy = k ■ r 
paraflel to the hne of sight and a 2-dimensional projection, k^ = k — fc||f , in the plane 
orthogonal to this line of sight. 

While the orientation of the Fourier mode with respect to the line of sight does 
not matter for the monopolc of the brightness = 0), it plays an important role for 
higher order multipolcs {£ > 1). Indeed their contribution to the radiation that wc 
receive depends on the orientation of k with respect to the direction of the flat sky. 

The flat sky approximation is independent from the resolution of the Boltzmann 
equation. We start by solving the first few moments of the emitting sources by 
integrating a truncated Boltzmann hierarchy, and only after this we use the line of 
sight method to obtain all the multipoles of the CMB. The flat sky approximation is 
only used in this final step and proves to be a very good approximation for £ > 10, since 
the error is much smaller than the cosmic variance, and reaches a percent accuracy 
beyond i = 100. Two subcases of the flat sky approximation exist [77]. The Limber 
approximation corresponds to cases where the sources vary slowly in time and over a 
wide range of distances along the line of sight. This case is not well suited for CMB 
except for effects like the late ISW. The second regime corresponds to the thin shell 
approximation. It relies on the hypothesis that the sources are located in a narrow 
range of distances. This is the one that we use in this article since it appears to be 
well suited for primary CMB calculations as illustrated in Fig. [TJ 

Let us also note that the procedure described in the previous section in order 
to consider a general orientation of k, ki and k2 with respect to the flat sky can be 
obtained as follows. First, we align the azimuthal direction k with the direction f 
of the flat sky which is chosen as the z direction of a Cartesian coordinate system. 
Two other axes x and y perpendicular z can then be defined. We choose ki and k2 
to lie in the plane spanned by z and x. From this particular configuration, we can 
reach a general configuration by first rotating around z by an angle a, then by rotating 
around y by an angle /3, and finally by rotating around z by an angle 7. These last two 
rotations enable to reach k from the direction z. In the notation Y(,'"(k) of Eq. (|4.6[) . 
we thus meant Yp{a,(3). 

4-4- Spectrum 

As usual, the temperature anisotropics can be decomposed in spherical harmonics as 



hence defining the coefficients The global isotropy of the Universe implies |78[[79] 
that the 2-point correlation function of Q can be decomposed as 



where n\n2i = cos 9 and the brackets refer to an ensemble average or to an average 
on the sky if Q is the observed temperature. The angular power spectrum can then 
be related to the ensemble average of the aim as 




(4.13) 




(4.14) 




(4.15) 



The CMB bispectrum from the cosmological perturbations 



15 



The angular power spectrum can be expressed in terms of the initial power spectrum 
and linear transfer functions [HI HH] • At lowest order in the flat-sky expansion, the 
angular power spectrum of the CMB anisotropies takes the form [TS] (see also |77| ) 

= ^/ drdr'dfcj|«;(i)(k,r)^*(i)(k,r')e-''=il('-'-')^^^:^^ (4.16) 

where r — rjo — rj, k± = {£ + l/2)/r and fc^ = fc^ + fc^, and where a star denotes the 
complex conjugate. If we neglect the integrated effects which occur at low z and affect 
only the largest scales, we can then approximate that all the signal in w(k, rj) arises 
from the last scattering surface whose thickness is defined from the visibility function 
g(r]) = — 'f'exp(— f). This visiblity function peaks at tlss which defines the center of 
the last scattering surface. In this thin shell approximation, we thus obtain that 

I / dr^W(k,r)e-''=il'-':^ , (4.17) 

now with k± = {£ + 1/2)/7'lss- Additionally, since modes of different m are in 
general not statistically correlatcc0, the phases in the spherical harmonics of (|4.6p 
are unimportant and we can take for practical purposes 



Cf « ^ ^ / dfc 



QTi^) ^ (-i)y ih/k) , (4.18) 

where P™ are the associated Legendre polynomials. 
4-5. Bispectrum 

In a similar way as the angular spectrum, one defines the bispectrum which, again 
under the assumption of isotropy, is related to the 3-point function of the aim by 

jeee / ^2 «3 



(a^imiO^^m^afgrnj) = Bg^g^}^ \ mi 7712 ms j ' (4.19) 

where the matrix denotes a Wigner-3j symbol and Bf^f^^ is the bispectrum. It has 
become more usual to use the reduced bispectrum bf^f^^ [111 ISO] which is defined from 
the bispectrum Bf^f^^ by 



eee ^fh h h \ (24 + 1)(2£3 + 1)(24 + l) ,eee .n^ 
In the flat-sky approximation [77], it takes the form 

^W3 = ^ J dfci||dfc2||dridr2dr3e-"^iii('^i-''3)e-''==^il(''=^^''^) (4.21) 

\ (i)n \ *(2)., , ■,P>S>{ki) P<S>{k2) , „ 
X w^'' {ki,ri)w^'' {k2,r2)w^Y i^i^^^^^s) + 2 sym. , 

where 

fcij. = £i/ri , k2i_ = £2/r2 , (4.22a) 

rir2 rir2 \ri J rir2 \r2 ) 

\ Note that in our case this is also the case because we consider scalar modes at first order in the 
perturbations. 
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Wc also remind that, since 1*3 = ki +k2; the parahel components are related according 
to A:3|[ = fcj^ii + fc2||- Similarly to the spectrum, this expression can be avantageously 
simplified in the case where the emitting sources peak in a narrow range of distances. 
More precisely, provided that the first order emitting sources are peaked in a narrow 
range of distance, which is the case if we neglect the late ISW effect, the following 
approximation holds out of the largest scales. In order to grasp how sharply peaked 
the first order emitting sources are, we plot them in Fig. [1] for £ = 0, 1, 2 for scalar 
perturbations (m = 0). In this thin shell approximation, the expression ()4.21|) can be 



Sao. S|o, Sio, k/t(|=10 




^00^ Suh ^Kir k//:(.c|-0.1 


' ' i - 


1.0 
0.8 






0.6 






0.4 






\ 

\ 


0.2 








0.0 










1 


I 


-0.2 








-0.4 





1 2 5 10 20 50 100 1 2 5 10 20 50 100 



Figure 1. First order scalar emitting sources 5^ witli ^ = 0, 1, 2 in respectively 
solid, dashed and dotted lines. The plot on the left corresponds to a mode which is 
well inside the Hubble radius around recombination time, whereas the plot on the 
right corresponds to a mode which is still super-Hubble around recombination 
time. In both cases the contributions are sharply peaked and this defines the 
last-scattering surface and validates the fiat-sky approximation. 



simplified to give 



Leee 



1 

2^ 



dfc]^l|dfc2 



P^{ki)P^{k2 



2 perm. (4.23) 



'LSS 



now with the simpler fiat-sky relations 

ki± = £i/r-Lss , k2± = e2/rhss , ^si = ia/r-LSS ■ (4.24) 

From the rotational invariance of the spectrum, we can also use the expression (|4.18p in 
the computation of w'^^^ and w^^^. Furthermore, given the definition of the multipoles 
with the factor iV^ |53j . the symmetries fci|| ~> — fci|| and fc2|| ~^ ~^2|| correspond to 
a — > TT — a, and since bf^^ is real, we can replace exp[ima(ki, k2)] by cos[ma(ki, k2)] 
in Eq. gUJ. 



5. Numerical implementation 

5.1. Width of correlations 

In the expression (|4.16p of the spectrum, the term multiplying e"'*^!! ^ is an almost 
constant in the regime fc|| ^ k±^ and is damped, at least by the power spectrum 

P<s> ( + when fc|| ^ k±. It follows that the integration over fc|| leads to a 
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window function [40j of argument (r — r') whose width is of the order l/k± ~ ''lss/^- 
A similar argument holds for the bispectrum and the integration over and fc2|| 
in Eq. (|4.2ip leads to two window functions, one for (ri — r^) with a typical width 
l/ki± ~ J'LSs/^ij and a second for (r2 — ^3) with a typical width l/fc2± — fLSs/-^2- 
We thus conclude that for a given mode £, the typical scale of correlation is given 
by tlss/^- This justifies why we can safely split the effects between primary and 
secondary effects, that is between early and late time effects. Indeed, at first order, 
the gravitational potential reaches a constant value until the Universe starts to be 
dominated by the cosmological constant. The integrated effects are then separated 
into an early effect, that is before the potential freezes, and a late effect when it starts 
evolving again. Since the correlations of physical effects must occur at approximately 
the same distance, then even if integrated effects occur everywhere at second order as 
the potential grows, they would only contribute in the CMB bispectrum by correlation 
with the first order early and late effects. As a result, when we are interested in the 
bispectrum, the second order effects can still be split into primary and secondary 
effects. In the numerical implementation of Eq. (|4.2ip . this implies that, as long as 
we neglect the late ISW, we only need to integrate over ra in a range of tlss plus or 
minus a few times tlss/^ to encompass the early ISW. 

5.2. Numerical integration 

The numerical integration can be sketched as follows. 

(I) We first integrate numerically the sources. Using the initial conditions defined in 
§ m we integrate numerically in Fourier space the coupled system constituted of 
(1) the Einstein equation, (2) the Boltzmann equation both up to second order 
in the perturbations and (3) the recombination equation up to first order in the 
perturbations. We first need to integrate the first order equations and then use 
their results to determine the sources of the second order integration. 

(II) We then compute the line of sight expressions (j4.17p and (|4.23p respectively for 
the spectrum and for the bispectrum, with their flat sky approximation. 

Let us emphasize the following details: 

• All our numerics assume the cosmological parameters as determined by WMAP- 
5 [60]. 

• The numerics requires to sample the Fourier space and for each Fourier 
configuration to integrate numerically in time. At first order, we only need to 
sample the norm k of the mode k, since we aligned the mode with the azimuthal 
direction. At second order, we also align k with the azimuthal direction. We only 
need to sample ki and k2 and the cosinus of the angle {fJ.12) between ki and k2. 

• The Boltzmann hierarchy has to be truncated in order to avoid spurious 
reflections. We use the closing relations defined in Ref. [7^ for both the first 
and the second order. The latter being defined only for scalar modes, we need to 
supplement them with the appropriate closing relations for the modes with m > 1 
that appear at the second order. These can be found in Ref. [80] and we report 
them as well in [Appendix D| 

• We use f max = 8 at first order and £max = 5 at second order for the truncation of 
the Boltzmann hierarchy. 



The CMB bispectrum from the cosmological perturbations 



18 



• The line of sight integrals for the spectrum and for the bispectrum are then 
sampled only up to 3 times the caracteristic scale of the correlation (i.e. ^lss/^) 
from the end of the last scattering surface. 

• We do not include the late ISW and our results hold only for £ > 10. 

• The code used runs in Mathematica and is freely available with documentation 
on the webpage [58]. The package CMBquickl contains all the functions used for 
the first order integrations including the transfer functions and the spectra, while 
the package CMBquick2 contains in the same way all the functions used in the 
second order integrations. 



5.3. Transfer functions 




Figure 2. Top left: and f'^) 

are in plotted in solid and dotted lines 
respectively. Top right: Vector perturbation B^^' (solid line); Bottom left: Tensor 
perturbation H^^^ (solid line). For all these panels, the asymptotic transfer 
functions in the limit a /acq ^ 1 are depicted in dashed line. Bottom right: Second 
order energy density contrasts for the radiation, Sj (solid), and baryons; Si, (long 
dashed). The tight coupling approximation which holds as long as t' /k ^ 1, is 
drawn in dotted line. 



Fig. [5] depicts the transfer functions at second order, using a definition similar 
to Eq. (|4.1ip . for the scalar, vector and tensor degrees of freedom of the metric, 
restricting to the particular configuration ki = k2 = &kcq and fii2 = —0.5. These 
transfer functions, computed numerically and without any approximation are perfectly 
compatible with the ones reported in Ref. [49] obtained from an approximation of 
the Boltzmann equation. This gives an a posteriori confirmation of the physical 
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understanding we have provided in Ref. [49] and we refer to this article for further 
discussion. 

Fig. [5] also compares the transfer functions to their asymptotic behaviours in the 
limit a/ooq 3> 1. These asymptotic behaviours correspond to the transfer functions 
computed in a purely matter dominated universe, and we have checked that they 
match with their previous estimations j40j . 

Note that the vector modes transfer function almost vanishes when ki — k2 for all 
Hi2 since in that case the only contribution arises from the coupling of the anisotropic 
stress contracted with gradients [SJ . Indeed in the harmonic decomposition of terms 
like XdiY, the terms for m = 1 and m = —1 vanish for symmetry reasons when 
fci = ^2. For such configurations, the contributions to vector production then arise 
solely from terms of the form X''^djY. As a consequence, we have plotted a Fourier 
configuration which escapes this restriction. 



5.4- Reduced bispectra 

5.4- 1- Relation between bispectra In § 14. 1[ we emphasized the different possible 
definitions of temperature. Each definition defines a bispectrum while it does not 
affect the spectrum. 

While the most straightforward would be to use the bispectrum of the temperature 
brightness Q defined in Eq. (|4.ip . the previous literature |3H] has focused on the 
bispectrum of the bolometric temperature T defined in Eq. (|4.2p . And, as we 
emphasized, to avoid the spectral distortion that affect the bolometric temperature, 
it may be more robust to use the occupation number temperature Tjy. Let us now 
discuss the relation between the bispectra of these 3 quantities. 

First, it is obvious from Eqs. (|4.ip and (|4.2p that 

bJZs = bZT" - 3 (^^1^^. + C^^Cis + Ce.Ce,) . (5.1) 
Then, we deduce from Eq. (|4.3p that the total second order information is carried by 
the combinattion of bj^^^f^"" and b^ef^J^ ■ It follows that 

■ TnTnTn _ jTTT _ t^TnTn _ .^TnTn _ t^TnTn /r rjN 

"111213 ~ "«l«2f3 "^1^2f3 %f3^1 "«3«1«2 ' ^'^■^) 

In conclusion, it is an easy translation to switch from one bispectrum to the other 
(once the evolution of y is known but this has been achieved in Ref. (S^). In order 
for our results to be easily compared to the previous literature, we will only consider 
the bolometric temperature. 



5.4-2. Examples of numerical results Representing the bispectrum completely would 
require a four dimensional plot. For the sake of clarity, we restrict our presentation 
to a few specific configurations. But once again any other configuration is accessible 
via the use of the packages made available in [58] . On Fig. |3l we show the explicit 
shape of the bispectrum that we have found and show the three contributions coming 
from the sources with m = 0, m = 1 and m = 2 respectiveljQ- We see that the scalar 

+ Note that the meaning of I for a source term is indeed ambiguous since for instance if we consider 
the gradient of a scalar, it contributes in the source with 1 = 1 and m = in an integrated effect. 
However if we perform an integration by part of the form n^djX = AX/drj — dX/dri then the total 
derivative contributes in the source term £ = m = in an effect on the last scattering surface, and 
the partial derivative contributes as well in £ = m = in an integrated effect. Actually, at first 
order, this is precisely what happens to the term n}di^'^^^ and that is why after an integration by 
parts it contributes to the Sachs- Wolfe effect and as an integrated effect. Furthermore, there is no 
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contribution, that is from sources with m = is dominant, especially on small scales. 
The top panel corresponds to equilateral configurations (with a varying size) and the 
bottom panels to squeezed cases (isosceles triangles) with length ratio of 10 for the 
middle panel and a fixed value h = 20 for the bottom. The curves are given as a 
function of ^3. 

6. Implications for the detection of the bispectrum 

The previous sections have allowed us to numerically compute different bispectra 
taking into account the full non-linear evolution of the cosmological perturbations 
and compare it to different types of primordial non-Gaussianities. Wc now turn, in 
§ 16.11 to the question of defining an equivalent /^^ then to the estimation of the 
signal to noise ratio in § 16.21 

6.1. Equivalent f^^^ 

Wc follow the extensive previous works [25l |35l [38l [39l |40l [71] in order to define what 
is the equivalent /^^j noted f-^-^, that the primary non-linear effects would induce. 

More precisely, let us consider an ideal experiment, that would be only cosmic 
variance limited, and a statistical estimator designed to detect the primordial non- 
Gaussianity of the type f^^ assuming a linear evolution of the cosmological fields. 

Then /^^ is precisely what such an estimator would actually measure because of the 
non-linear dynamics. In summary, for a given choice of model for the primordial 
couplings, is such that 

/ fL-l. \ I fl" = ^ \ 

Signal with = Signal with . (6.1) 

\ linear evolution / \ non-linear evolution J 

According to this point of view, the signal coming from the second order evolution 
represents a spurious signal, a noise, and /^^ is the amplitude of that noise. Of course 
from an observational point of view, this is not a genuine noise and is not related to 
any experimental imperfections. It may also have some interests on its own! In the 
future, it will have to be incorporated consistently in the design of the estimators of 
primordial non-Gaussianity. 

Being mostly interested in the equilateral and the local configurations, we report 
in Fig. [4] /j^j^ for these two configurations as a function of the maximum multipole 
^max used. 

We have also plotted on Fig. |4] the f^^ arising only from the scalar (m = 0) 
purely second order source. As expected from our previous analysis [49] and from the 
numerical results of the previous section, these are clearly the dominant contributions. 
When comparing with the f-^-^ coming from the secondary non-linear effects as reported 
in Ref. [39] , we observe that the primary effects are of slighlty less importance for the 
local case while they are far more important for the equilateral case. For the local 

physical meaning in splitting the contributions of the sources into purely second order terms and 
quadratic terms as in Ref. 1381 since this decomposition is not gauge invariant. Additionally, there 
can be large contributions of opposite sign which thus nearly cancel as for instance on large scales 
(see S l7.3l or Refs. [321 140j ) so that picking up a single contribution may be misleading and give an 
overestimation of the signal. With the same argument, the coupling of the first order collision term 
with the perturbed ionization fraction which has been studied in Refs. [361 1371 [82] is also not a gauge 
invariant quantity and we do not try to reproduce these results for this reason. 
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Figure 3. Comparison of the bispectrum induced by the non-linear dynamics 
(thick solid line) to a primordial bispectrum with = 5 (thick dashed line 
for local type, thick dotted line for the equilateral type). We also plot the 
approximation of Ref. I49| that is considering only purely second order scalar 
sources in thin solid line, and the contribution from sources for m = 0,1,2 
in respectively thin dashed, dotted, and dot-dashed lines. We considered 
three different configurations in I space, from top to bottom: an equilateral 
configuration with £1=^2= ^3 and then two squeezed configurations with 
10£i = £2 = ^3 and = 20 with £2 = h- 
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Figure 4. Comparison of /^^ (solid line) induced by the non-linear dynamics 
to /fjL obtained by including only the scalar (m = 0) modes (dashed line) of 
the purely second order sources, which corresponds to the approximation of our 
previous work 1491 . We also plot /^j^ (dot-dashed lined) induced by lensing-ISW 
secondary effect. Right panel corresponds to local type couplings and the left 
panel to equilateral type couplings and the icr detection limit of f^^ is depicted 
in dotted lines for the local type couplings case. 



case, the contributions to /^j^ from primary and secondary effects have the same sign. 
The sum of the primary and secondary effects is then beyond the 2a hmit of detection 
when ^inax > 1500. This is certainly above the detection limit of Planck since it 
uses multipoles up to ^max — 2000 and it should thus be taken into account. The 
primary effects alone would bias the measurement of the local type non-Gaussianity 
by A/* ^ 5. 



6.2. Signal to noise 

We can now report the signal to noise ratio (S/N) as a function of the maximum 
multipole iraax uscd in the estimator described in Refs. [321 Hi]- Fig. [5] represents the 
signal to noise ratio of the primordial bispectra, both for the local and equilateral 
cases, as well as the signal to noise ratio of the bispectrum generated by non-linear 
effects. We also superimpose the signal to noise ratio of the bispectrum generated by 
ISW-lensing correlations (see § 17. 2[) . We found that the signal to noise ratio of the 
bispectrum generated by the non-linear dynamics reaches unity for ^max = 2000. 



7. Analytical understanding of the limiting cases 

The aim of the section is to derive analytical formulae, a priori valid in asymptotic 
regimes, that offer theoretical insights on the underlying physics. We first ignore the 
contribution of neutrinos, and thus set = 1, since they play only a marginal role it 
will greatly simplify the expressions obtained. 

After rewriting the Boltzmann equation in § 17.11 we first focus on secondary 
anisotropics in ? I7.2I We then consider the large angular scale limit in ? I7.3I (considering 
separately local and equilateral configurations), the small angular scale limit in ? I7.4[ 
and finally a mixing between small and large scales in § 17.51 
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Figure 5. The signal to noise ratio as a function of the maximum multipole 
^max for an ideal experiment. The signal to noise ratio of the total bispectrum 
generated by non-hnear effects (sohd hne), of the local type primordial bispectrum 
(dashed line) and of the equilateral type (dotted line) when f^^ = 1. We also 
plot the signal to noise ratio of the bispectrum due to the ISW-lensing correlation 
in dot-dashed line. 




7.1. Rewriting the Boltzniann equation non perturbatively 

Since the Boltzmann equation is solved numerically by using the line of sight approach, 
it is convenient to rewrite it in terms of the total time derivative along geodesies. It 
has been noted that defining 

exp[2$BRM] = 1 + 2$ , exp[-2«'BRM] = 1 - 2* , (7.1) 
simplifies its expression |85| . In particular, the Boltzmann equation, written non 
perturbatively (that is in such a way that it is valid up to any order in perturbations), 
can be recast as 



djy (771* arj h dyy 

In this expression, E stands for the scale factor multiplied by the energjQ of a photorQ, 

I(x, 77, n') is the reduced brightness of the radiation, that is the brightness divided 

- -0 

by its average value (this implies that I = Iq = 1). Note also that d/d?7 is the total 
derivative on the perturbed space-time which includes the convective derivative, but 
not lensing, i.e. defined on the perturbed space-time by 

d _ d dx^ d 

d?7 dr] di] dx^ 

While being very general, E cannot be specified so easily in a non pcrturbative manner, 
and neither can we obtain a non pcrturbative expression of the collision term. However, 
up to second order we obtain 
1 AE _ 
E'd^ ^ 



(7.4) 



* In Ref. 1531 E is noted ap° and in Ref. 1351 it is noted ap. Since this is the energy of photons, we 
use this more obvious notation. 

(j See Ref. |53| for details about the observer which actually measures this energy. 
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where we must bear in mind that the vector and tensor modes are only second order 
quantities. As for the colHsion term, it is reported in |Appendix A.2] up to second order 
as well. 

However, we can go a few steps further in our non-pcrturbative reformulation 
of the Boltzmann equation. Indeed, we can integrate by parts Id(biE)/dr], and 
use the Boltzmann equation itself to replace (In _B)dZ/d77. We then use that 
(ln£')d(lni?)/d7/ = l/2d(ln£')^/d?7 and iterate this method and finally resum all the 
terms to obtain 



d 



. di dn' i 



d?^ (?n* d?7 

If we expand the photon geodesic equation in perturbations according to 



C[I]+f'I . (7.5) 



dx^ dx^ ^ / dx^ " 
drj dri \ dr] 
we can expand accordingly the total derivative as 



(7.6) 



— - — ^ — = — + —— (7 7) 

d?y drj \d7] J ' dr] dr] drj dx^ 



An integral solution of Eq. (|7.4|) can be found immediatly by integrating the 
Boltzmann equation over the background geodesies. We obtain finally that the 
observed brightness today in the direction of observation — is 

E-\r],)T[,„-n^)^ j dr]'-gir]')E-^[c[lf—^+iyj e'^E'^^^ 

This expression is valid up to any order in the perturbations, but the difficulties are 
hidden in the computation of E and in the expression of the collision term C \L\ . The 
first integral on the right hand side encodes the intrinsic anisotropy. As expected, 
it is weighted by the visibility function and adds up the contributions from the last 
scattering surface. The second integral contains the lensing effect while the last one 
is the time delay. This latter contribution has already been computed [53] for the 
spectrum and for the bispectrum, and its effect is by far negligible for geometric 
reasons. However, the lensing effect is both a primary effect and a secondary effect 
since it contributes from the last scattering surface up to now. Since dn^ /dr] vanishes 
at the background level, we can use the general integral equation (j7.8|) to express it 
in terms of the intrinsic anisotropics by iteration. We detail in the next section how 
this gives rise to a bispectrum when coupled to secondary linear effects. 



7.2. Secondary effects 

After having reached a constant value in the matter dominated era, the first 
order potentials remain constant until the cosmological constant eventually starts to 
dominate the matter content of the Universe. This gives rise to a late integrated effect 
(which enters in E) known as late Integrated Sachs- Wolfe effect (late ISW). Given 
the arguments of § 15. li on the width of the correlation in the bispectrum, the intrinsic 
second order anisotropy will not correlate to the late ISW. However the lensing term 
can correlate to the late ISW since it is also an integrated effect. The lensing term is at 
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least a second order term, and if we make use of the first order Boltzmann equation to 
replace for dl^^^ /dn\ it can be recast in (see next section as well as Refs. [84 1 186 1 187 ] ) 

V">(iV„.i(7?o,-riO, (7.9) 
where 0*^^^ is the lensing potential whose expression we give further in Eq. (|7.18p . 
Note that in order to obtain this expression, we have to neglect the lensing inside the 
last scattering surface as well as the lensing of the late ISW. Both approximations are 
justified from the geometrical properties of the lensing effect since it vanishes when the 
lens is either at the source or at the observer position. Under these approximations, 
the lensing term is separable into a lensing term 0^^^ and a lensed term The 
computation of separable terms is straightforward since the lensing potential correlates 
with the late ISW effect and the lensed source correlates with the temperature emitted 
on the last scattering surface. Strictly speaking the bispectrum arising from the 
lensing term is thus a combination of primary and secondary effect Jpjl . The reduced 
bispectrum obtained is given for instance in Refs. [221 [Ml Ell [EH UHl [HZl [HH] and reads 

Ki.i, = \ Mi2 + 1) + 4(^3 + 1) - h{h + 1)] Cf^'^-'CI'' + 5 sym., (7.10) 

where the prefactor comes from the angular gradients on (f)^^'^ and This method 
to derive the correlation of the lensing term with first-order secondary effects is 
general [551122]. Reionization for instance will create a late Doppler effect, and the 
effect of the non-linear structures formation (Rces-Sciama, Sunyaev-Zel'dovich...) can 
also be mapped to additional secondary linear effects (in an effective manner) in order 
to apply this method [41] |89] and compute their bispectrum. 

7.3. Large scales 

7.3.1. Generalities On large scales, that is for modes such that krji^gg <C 1 where 
Vi^ss is the average time of recombination, we want to derive the equivalent of the 
formula O'-"^' = ^$lss that allows to explain the Sachs- Wolfe plateau of the angular 
power spectrum. This was already obtained by following the photons geodesies in 

Refs. [saissiiio]. 

We assume that the last scattering surface is far in the matter dominated era, 
that is long after the matter-radiation equivalence. We also ignore the effect of the 
cosmological constant. The Bardeen potentials and have then reached a 
constant value that we label by LSS. On large scales, this is even true at any order 
in perturbations. The width of the visibility function can also be neglected on large 
scales, and we approximate g{ri) by a Dirac distribution centered around ?7lsSi and 
e""^ by a step function which is unity for rj > t^lss. We can also neglect the velocity 
fields, that is the velocity moments Vm of baryons and cold dark matter, and the first 
moment of radiation and neutrinos, that is Z™ and A/"™ since they are suppressed 
by fcr/Lss. Higher order moments can also be neglected since they are suppressed by 
higher powers of fcyyLSS. 

A first simplification in Eq. ()7.8|) lies in the fact that the collision term has no 
monopole on large scales. We further assume that we can neglect the vector and 
tensor modes and the integrated terms in the derivative of E given by Eq. (|7.4p . In 
that simplified case we have 

E = e"*^™ . (7.11) 

ft Pure secondary effects appear only at large scale and are thus less interesting since the cosmic 
variance limits the information extracted from the observations. 
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The general solution obtained in Eq. (|7.8p then reduces to 

rno a-f j„i 



17LSS 



dn^ dr] 



<5 pi^ 3/ Ie4*«™)d77', (7.12) 



drj 



which is also valid non perturbatively (that is up to any order in perturbations) 
provided that the approximations hold. We will discard from now on the last integral 
coming from the time-delay as it is negligible [53]. Restricting to an expression valid 
only up to second order, we obtain (restoring the potential $) 

I(r/o,-n') -4$lss+2:olss+4(2:lss + $lss)$lss- / TTT-^dr;'. (7.13) 

Note that in the previous expression, following standard practice, we have omitted the 
terms which depend on the observer space-time location as they are not observable. 
At first order, we obtain the conservation relation 

i^5!!j±f^.*.(.,-«s)(e"> + *'«), (7^14) 

S]j is the one-c 

result 



where Sj-, is the one-dimensional Dirac distribution. We thus recover the textbook 



e(i)(77o, -n') = [e'J'^ + (r^Lss) , (7.15) 

where we recall that G is defined in Eq. (|4.ip . We also recall that the right hand side 
is evaluated at t/lss on background geodesies defined by the direction and that we 
have dropped the potential evaluated today since it is not observable. 
At second order, we have 



e(2)(r;o,-n') = e°(Lls + <sU2 



no 50(1) /^ni 



-"lSS + ^LSS 

(1) 



The first line of the right hand side represents the intrinsic anisotropy while the last 
line contains the lensing effect. Using Eq. ()7.14p and the fact that the gravitational 
potential does not depend on the direction of propagation but only on the space-time 
coordinates, the lensing term can be recast as 

-/::^(w)'"^''-"^*'"'"'Hs...ss), 

where we have used the first order expression for dn^/dr] (given in our notation in 
Ref. [53]). and where we have defined the lensing potential by 

. - r [*"'w)+i-'w)]w-.>Lss) 

the geometic factor being simply the usual combination of angular distances in a 
spatially Euclidean universe. 
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For adiabatic initial conditions, and assuming that the Universe is completely 
matter dominated during the recombination, we obtain the following initial conditions 
for the temperature anisotropics 

which can be obtained easily [40l [55] using Ilss = exp (— 8$B|uy[/3) . Furthermore, 
we can use that in the matter dominated era the second order gravitational potential 
is given on large scales by 

^(2) ^ 2$(i)2 + 2A-i(9/$(i)9^4>(i)) - 6A-2a^9-'(aj$(i)a,/$(^)) , (7.20) 

when = — 1 . This enables us to obtain a complete expression of the second order 
temperature anisotropics on large scales in the case where the approximations made 
hold. 



7.3.2. Local configuration In the local configuration, the primary non-Gaussianity 
contributes mostly to squeezed configurations of the bispectrum. that is when one of 
the modes is much smaller than the other two, for instance if fci <C A:2. We are thus 
interested in the contribution of evolution in that same configuration in the £ space in 
order to assess the contribution of non-linear effects to the total CMB bispectrum. The 
contributions of the purely second order vector and tensor modes are suppressed by 
ki/k2 and can thus be ignored. The same happens for the integrated effects involving 
$'2^ and ^''2'' . The equation (|7.16p is thus valid in that case. Furthermore, the 
lensing term [40l [90] can also be neglected in this large scale limit, and we are thus 
left with the intrinsic anisotropy. For the local configuration we read from (|7.20p that 
$(2) = 2$(i)2^ g^jjj ^j^^g obtain 

9(1) ^ l$(i) ^ e(2) = l$(i)2. (7.21) 
o y 

We recover the results of Ref. [ID] for the bolometric temperature T, since the 
expression obtained reads precisely up to second order 

(T/T)* = I = (1 + 46) = cxp[4$(i' /3] . (7.22) 

Should we need to compute the temperature with a non- vanishing f^^'^' , that is 
with 7^ —1, the previous result would be modified by replacement of by 
- + 1) ($(1)2 _ ($(1)2^^ and the arguments exposed in Ref. [iO] for the 
interpretation of the result would remain unchanged. Note also that this ideal regime 
is never reached since the last scattering surface is not deep in the matter dominated 
era. Since the first and second order potentials undergo a change of value between 
the radiation dominated era and the matter dominated era, the relation $*^2) = 2$^^^2 
is violated around equivalence and in particular on the last scattering surface. This 
fact implies that the limit in which Eq. (|7.22p was obtained cannot be seen from our 
numerical results. 
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7.3.3. Equilateral configuration In the equilateral case, the previous approximations 
do not hold. Indeed the lensing term and the integrated contributions of the 
scalar, vector and tensor modes cannot be ignored since terms such as A~^d^^di^ 
are negligible in the squeezed limit but not in an equilateral configuration. Since 
these effects do not correlate only at the last scattering surface, a proper numerical 
integration of the bispectrum expression (|4.23|) has to be performed. They are 
computed in details in Ref. [40] . 

7.4. Small scales: Equilateral configuration 

On small scales, another approximation scheme can be used in order to assess the 
order of magnitude of the evolutionary bispectrum. This has been presented in detail 
in our previous analysis j49j , the validity of which has been confirmed by the numerics 
of the present work. It relies on the assumption that all modes are on sub-Hubble 
scales, and it thus holds only for non-squeezed configurations, that is for configurations 
of the equilateral type. To summarize, it is based on the fact that any perturbation 
mode of the cold dark matter starts collapsing as soon as it becomes sub-Hubble. The 
non-linear gravitational potential is thus predominently shaped by cold dark matter 
component. The photon-baryon plasma then develops acoustic oscillations, that are 
imprinted on the last scattering surface, with a forcing term given by the non-linear 
gravitational potential since on very small scales the purely second order Sachs- Wolfe 
effect is driven towards 



and other effects can be neglected. The convergence of this approximation onto the 
complete computation is shown on Fig. |4l and we see that it is responsible for most 
of signal on small scales. In particular, we conclude that the approximation is good 
as soon as ^ > 500. Note also that it is in very good agreement with the full analytic 
estimation of Ref. [35| which also find /„l — 5 for equilateral type non-Gaussianity 
when ^max = 2000. 

7.5. Mixing large and small scales 

The fact that the convolution (|2.8p performed on the source terms is invariant under 
the exchange fci -(->■ k2 makes it very difficult to find analytic solutions to the evolution 
of perturbations in that case. To see this, let us consider a contribution to the 
second order potential which would be quadratic in first order variables, that 
is of the form [Xy](k) [using the notation of (|2.7p ]. The first order quantities X 
and Y can be related to the primordial first-order potential through a transfer 
function, these functions being defined for instance by X{ki) = Tx(fci)$^^''(fci) and 
Y{k2) = '?V(fc2)'i''"^'(^2)- The contribution to the bispectrum in $ of such quadratic 
term is then found to be 



([Xy](k)$(ki)$(k2)) = 4(k+ki+k2) [Tx(fcl)Ty(fc2) + fci O fc2] P$(fci)P*(fc2) .(7.24) 



We thus notice that if we want to find an approximation in a configuration ki <C /s2, 
we have to consider an approximation in that regime for both Tx{ki)TY{k2) and 
Tx{k2)TY{ki) . For such configurations, there is thus no simple approximation scheme 
that could be used and one has to rely entirely on the numerical results. This is 
typically an issue in order to estimate the non-Gaussian signal for local configurations. 
Indeed most of the signal in the convolution comes from configurations where one 



(7.23) 
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of the mode is much smaller than the other (for instance fci ^ fc2). And for this 
type of configuration, k2 might corresponds to sub-Hubble scales if k2 » fccqj but 
ki must also be much smaller and thus super-Hubble. One can however estimate 
the non-Gaussian signal for local configuration by comparing numerically the purely 

second order Sachs- Wolfe effect i (^6q^^' +$(2'^ (ki,k2) with fci < k2, to the first 

order one ^Oq'^'' + (fc), where we recall that k = ki + k2. For two types of 

squeezed configurations where the large mode ki = 15fcoq corresponds approximately 
to £ ~ 2000, we plot these two effects on Fig. [6l In order to stress that the purely 
second order effect is larger than the first order one, we multiply the first order 
effect by a factor 5. We note that the first and second order Sachs-Wolfe effects 
are oscillating nearly in phase with opposite sign. This result suggests that, in such 
squeezed configurations, most of the non-Gaussian signal is carried by Qq^'^\ And, 
given the sign definition of in Eq. p.22|) . the second-order acoustic oscillations 
have an amplitude such that it leads to /^^ — 5- 



1/2(^4+*)'-' t, = 15fa=15i„| /ii2=0.5, 5x(V4+4>)"'(k) l/2{S^/4+>S>P k, = l5k2=l5k,^ /Ji2=-0.5, 5x(V4+*)"'(k) 
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Figure 6. In continous line wc depict ^ (^©o'^' -|-<3>^^'^ (ki,k2) and in dashed 

line 5 X (Bq'^' -|- (fc). The vertical shaded area corresponds to the last- 

scattering surface, that is to the period during which the vast majority of photons 
are emitted (99%), except for those emitted during the reionized era. The ratios 
between the modes fci, fc2 and fceq is kept fixed and the angle ^12 takes the values 
0.5, —0.5 on respectively the left and right plots. We notice that the purely second 
order Sachs-Wolfe effect is nearly in phase with opposite sign with respect to the 
first order one during the last-scattering surface, and approximately 5 times larger 
in magnitude. 



8. Conclusion 

This article presents a complete investigation of the imprint of the non-linear dynamics 
on the CMB bispectrum. The calculations were carried with all the matter fields 
of the standard cosmological ACDM model included with three families of massless 
neutrinos. The numerical calculations make use of a full numerical integration 
of the coupled system of the second order Boltzmann and Einstein equations for 
both the photons, with their polarization, and the neutrinos. Furthermore, line 
of sight integrations include first order effects in the recombination history. The 
initial conditions correspond to adiabatic initial conditions with a vanishing intrinsic 
primordial non-Gaussianity (in the context of standard single-field infiation, that 
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implies that a contribution of the order of the slow-roh parameters has been neglected) . 

The numerical integrations were done with the cosmological parameters of the 
best fit model derived from the WMAP data [50]. This article is focused on the 
bispectrum of the temperature anisotropics. We have been forced however to define 
the temperature we use, namely in this article the bolometric temperature, since 
second order effects are bound to induce spectral distortions. This effect has been 
described in details in Ref. [57], but we do not expect though that a change of 
definition for the temperature, like the occupation number temperature, would change 
significantly the conclusions we have reached. This work also demonstrates that the 
second order Boltzmann equation can be exactly integrated numerically and be used 
to produce bispectra. The resulting shape and amplitude of those bispectra is the 
result of intricate phenomena. It is possible though to obtain theoretical insights into 
peculiar cases, at small or large scale for instance. Those results confirm in particular 
that, at small scales, the major mechanism at play is the impact of the gravitational 
coupling of the dark matter potential during the matter dominated era as it had been 
put forward in Ref. [49] . 

It is obviously difficult to grasp those results in details. In order to be able 
to compare the amplitude of those effects to primordial couplings, we have defined 
and computed equivalent f,^^ parameters, f,^^. They are defined in such a way that 
it is the signal a statistical indicator designed to measure primordial /^^ would get 
from the amplitude of the temperature bispectrum. We have found that for both 
primordial non-Gaussianity of local or equilateral types we have f^^ ~ 5 for £ above 
500. When compared to secondary effects, namely ISW-lensing couplings, primary 
effects are found to be of comparable amplitude. The former are however more efficient 
in producing a /^^ the local type ; the situation is reverse for the equilateral type 
for which the signal comes predominantly from the primary effects. 

Evaluations of signal to noise ratio however show that the non-Gaussianity 
induced by the primary second order effects in the temperature field alone can only 
be marginally detected by the Planck mission. Secondary effects are more likely to 
be detected. This is however the first ever explicit and complete computation of 
these effects and even though the concordant model does not offer a good chance 
of detection, it might be a good way to put constraints on alternative cosmological 
models. 

We remind that the numerical tools used in this article are freely available and 
can be downloaded at [58j so that many other bispectrum configurations and transfer 
functions can be investigated at will. The code makes uses of a flat sky approximation. 
It is accurate enough for our purpose and this approximation does not interfere with 
the resolution of the Boltzmann equation. The code should also be complemented 
with secondary non-linear effects, that have been investigated in Refs. [39 ] [4T ] [88] . 
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Appendix A. Governing equations 

Appendix A.l. Einstein equations 

We use the standard definition that for a species labelled by a, fia = Pa/{3H^). When 
unspecified, a runs on r, 6, c that is on photons, neutrinos, baryons and cold dark 
matter. We also use the definition k = 8ttG. In the expressions of the quadratic 
sources, it is implicitly meant that the two perturbation variables in each quadratic 
term are first order variables. We thus omit the order superscript in that case to 
alleviate the notation. We have a set of four scalar equations (see Refs. pSl [5^ 

A*(2) _ 3^^(2)' _ 3^2^(2) - I ^ p^2) ^ ^ 

a 

3 



S3 , 



^-(2) _$(2) 

a 

Their quadratic source terms read respectively 

Si = -g^-A^-- 397*9^* -3*'2 + K^p,(l + u.,)9jKaVa, 

a 



($ _ 5»)A$ + 8-H($ - + 4($ - , 



4*2 _ 2'I>2 



A 



+ ^)d'{^ + $) + 2*A$ + J2 ^PaC^ + Wa)dlVad'Va 

+ 3(AA)-i9/aj + + $) + 2*9/9,7$ 

+K^p,(l + i(;„)9^Ka'Va 

a 

5*4 = 2H$^ - 4*^-' + 2dY\^'di'^) 

+ Kj2Pad7' [(1 + U.a)*a/K - (1 + cl^SadlVa] . 
a 

These expressions arc easily computed in Fourier space using that [53] 
[9^y(i)a^y(i)](k) = -lc{ki.k2V^^\ki)v'-^\k2)} 

ln=-l 



(A.2) 

(A3) 
(A.4) 



(A5) 



(A6) 



(A7) 

(A.8) 

(A.9) 
(A.IO) 
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= -/C{ki.k2i;^^'(fci)«(^'(fc2)} (A.ll) 



since 



in function of the multipolar components of the modes fcj"'' and fcj , we can use the 



[a/y(i)](k) = ifc4^'(A;), (A.12) 

^;W(k) = -fc('")^;W(fc), (A.13) 
1 

ki.k2 = kik2i - (-1)"^^"^^^"^ > (A.14) 

n=-l 

where fc/ = kj/k and fc^ = kjk^ . The notation F denotes here the scalar part of the 
velocity field in the coordinate frame, whereas the moments Vm are taken in a local 
Minkowski frame. See Ref. [53] for the precise definitions of these decompositions. 
Note also that in order to express terms like 

3{AA)-^didj - A-^di^d'<^ , (A.15) 

le multipolar components o 
identity (see [Appendix A^ for definitions) 
1 

J2 o4^2 fci"^fc^"^ = 3(k.k0(k.k2) - ki.k2 . (A.16) 

n=-l 

As for the second order metric vector modes, we can determine them from the 
constraint (that we only report here in Fourier space) 

- yS^?) =Kj2pa{l + Wa)v^S + (A.17) 

a 

with 

Sy=IC |4*(i)'(fci)4™'$^'nfc2) (A.18) 

+2.^p41 + wa) - <i>(^)(fci) - <i^(fc2)| . 

Finally the second order tensor modes are determined by a second order 
differential equation which describes how they are sourced by quadratic terms, 

i7(?)"(fc) + 2HH(?)'(fc) + k'H^^ = I + -P.^f^^'^) (A.19) 

+ liKl^lC |fc^^'fcW[<i>(l)(fci)$(^'(fc2) + Vl/(l)(fci)*^^'(fc2)] + E -pA'}ihWa}{k2) 

Appendix A. 2. Boltzmann equation for radiation 

Appendix A. 2.1. General remarks As mcntionncd in 12.31 (see also Ref. |53j ) it is 
sufficient to specify the purely second order part and the quadratic part of an equation 
to fully specify it at first and second order. The Boltzmann equation, which is formally 
written 

C[X]=C[X], (A.20) 
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where X stands for X, £ and and where C and C are respectively the Liouville 
operator and the coUision operator. We decompose the operators in the form 

C[X, A'(i), A-f^)] = C[X] + £(1) [X, A'(i)] + i (£^2) ^(2)] _^ ^(i)(i) -^(i)-^ ^ ^p^ ^l) 

and similarly for C. The details of the derivation with our present notation can be 
found in Ref. [53] (see also Refs. [851 |92])- Contrary to what has been performed 
in Ref. |53j we choose to report here the Boltzmann equation in the form £[] = C[] 
rather than in the form £*[] = C'^\\^ in order to facilitate the comparison with existing 
literature. We will use the following definitions that will simplify the notation 



/(£2_^2)(£2_g2) 



P 



±1 



2^2 



±1 



/(£±m)(£±TO-l)(^2_g2) 



2£2 



xrn 

Ap — 



1 /(£+ 1 ±m)(£=F"^) 



(A.22) 
(A.23) 
(A.24) 
(A.25) 



We report first the link between the moments of radiation and the fiuid description 
involving the energy density and the velocity of radiation. The link is similar for 
neutrinos. More details can be found in Refs. [53l |67]. Note that we define the 
multipoles I™ from X/I so that it is dimensionless, and we do the same for the 
polarization and neutrinos multipoles. At first order we have 



Ml) 



(fc) 



and at second order we obtain 



r0(2) 



(k) 



j-r,i{2) 



(k) 



Pi- 

^E(k) 



-2/c|- ^ (-l)"<"(klR-n(k2) 

Note also that since w^(2) ^ ^/(2) _ 2*(i'F^(i) for all species g3], then 



v^i\k) 



vi'\k) = o, 



v^^\k) = - fcT/(2)(fc) - 2/C {¥^\ki)k^"W^^\k2)} 

V^^\k) = (fc) - 2/C {*(lHfcl)fc2*V(l)(fc2)} . 



(A.26) 

(A.27) 
(A.28) 

(A.29) 
(A.30) 
(A.31) 

(A.32) 



Appendix A. 2. 2. First order At first order we choose to align the Fourier mode 
considered k with the direction with respect to which the moments are taken, and the 
dependence in k becomes only a dependence in its magnitude k. The set of equations 
obtained at first order is (dropping the obvious dependence of all quantities in rf) 



£(i)[i]?(fc) =i;"(fc) + fc 



2^ + 3 



'^'n^,{k)-^iU{k) 



2i~ 1 



S'^m'ik) - 6j4k<P{k) 



(A.33) 
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£W[£]°(fc)=£;°(fc) + fc 



2e + 3 



21-1 



C(i)[I]°(fc) = f \ -Il{k) + 5°X{k) + 4<5i«o(fc) + 5j- [l°(fc) - ^/6f °(fc) ^ , (A.35) 



C(i)[£]r(fc) 



Z2"(fc)-\/6£2°(fc) 



(A.34) 



(A.36) 



The first order magnetic modes are not excited since the first order vector and tensor 
modes are ncghgible, so wc did not report them in the above equations. Additionahy, 
the intensity and electric multipoles are only excited for m = 0, that is the reason 
why we also only reported this case. 



Appendix A. 2. 3. Second order At second order wc can only align k with the 
azimuthal direction of the moments, but not ki or k2 at the same time. However 
for a first order quantity the components for a mode in a given direction ki can 
be obtained by rotating the components X° obtained when we had decided to align 
this direction with 63. Namely this rotation leads for a mode k to 



An 



y*^™(k)xo (fc) , 



(A.37) 



21 +V 

and in particular for the first order velocity 

z;„(k) --fc(")«o(fc). (A.38) 
We finally obtain for the linear terms in the second order Liouvillc operator 

r O^Ttrrn I 7('n 
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(A.41) 
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As for the quadratic terms, wc obtain 

[i]r(fc) = 41^!,,,,,[i]r(fc)+4e^li4Hr(fc) 

n— — 1 
1 71 I T^- rn 

+ E g^ir-r(k2)4fci"^$(fci)-4vi/'(fcoxr(k2) 

n=-l 

+4Sl [^>{ki) - $(fci)] $(fc2)fc^"^ - 8(5,°*(fci)*'(fc2)} , (A.42) 

with 

<T-,.UymT{k)^ -2/C ^ <^l™+"(k,)fc(-")[$(fc,) + vl;(fci)] 

71— — 1 

+ 2^ E f^^r-T"(k2)4"' I^flfci) + , (A.43) 



n=-l 



+ 2^ E f^^rr (k2)(^ - mi) + ^'(fci)] . (a.44) 

n— — 1 

For the electric-type polarization, wc can also separate explicitcly the lensing and 
time-delay contributions 

+ - E #^^0"(k2)4fcr^$(fci) 
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n 1 1'™ 
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with 
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+ 2^ E |^^rT"(k2)4"' + , (A.46) 
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2^ E |^^rr(k2)(^ - l)fci") [$(fci) + *(fci)] . (A.47) 
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For the magnetic-type polarization, we obtain (note that we have corrected for 
mistakes in Ref. [53] which were pointed out in Ref. [55] ) 



^ ^ 71 — ~1 

1 



^ ^ n— — 1 



(A.48) 
(A.49) 



-2/C 



(^+1) ^ 



(A.50) 



Note that the lensing terms clearly differs from Eq. (44) of Ref. [38] , and this affects 
the results presented in Eq. (52) of Ref. [3^ (and subsequently this error spreads 
as well into Eq. (450) of Ref. [IS]). Indeed the harmonic expansion of the lensing 
term (|7.9I) should be treated as in section IV. A of Ref. [Mj. There, the screen projector 
which is implied in our notation V" (see Ref. [53] for details) is not split 



before performing the harmonic decomposition as done in Ref. 



into (5j and n 

but is instead used to ensure that the object manipulated in the lensing term remains 
tangent to a two-dimensional sphere whose radius is parameterized by n'. This enables 
to work with spinned spherical harmonics as in Ref. [84] and to derive the result almost 
immediately. Additionally, the multipoles of the radiation brightness {Ag('k) in the 
notation of Ref. [S^) are not rotated correctly when the azimuthal direction, with 
respect to which the spherical harmonics are defined, is shifted from the axis ki or k2 
of the first order calculation to the axis k of the second order calculation. Indeed, this 
change of reference axis leads to a transformation given by Eq. (|A.37p . Though this 
is taken into account correctly for the baryons velocity (for instance in the eleventh 
line of Eq. (44) of Ref. [55]), it is clearly forgotten for the brightness moments in the 
last term of the fourth line or in the last line. These two reasons explain why our 
treatment of quadratic terms is substantially different from Ref. |38| . 
Using a similar method applied to the collision term leads to 

r0(2)/ 



c(2)[i]r(fc) = 



(A.51) 



where P"<'^\k) is non- vanishing only if — 2 < m < 2 and is defined in that case by 



10 



(A.52) 



For the electric and magnetic type collision terms, their second order linear 
components read 



(A.53) 
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The quadratic terms arc then given by 

cw(i)[i(i)]r(fc) = 

r 1 „^^-m 



(A.55) 
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Appendix A. 3. Fluid equations 

For baryons and cold dark matter we only need the two first moments, that is the 
continuity and Euler equations, since there is no anisotropic stress. The continuity 
equations are 
-, / 

+ fcv^^^ (fc) - 3*'(i) (fc) , (A.58) 
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The Euler equation reads 
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In the case of cold dark matter, the right hand side of all these equations vanishes 
since it is collisionless. 



Appendix A. 4- Perturbed recombination 

The fraction of free electron at the background level is obtained from a Boltzmann 
equation where the collision term encodes all the relevant processes in the 
recombination. We shall not detail any of it here and it can be found in 
Refs. [ini [701 |93l |94] (see also the different method of Ref. [95|). Since the interaction 
rate is given by r' = aXcTT-ofT, where rig is the number density of electrons (free and 
non-free), (Tt the Thomson cross-section, and Xq the fraction of free electrons, its 
perturbation is given by 



' (1) ^(1)^ 



(A.62) 



Since electrons and baryons are tightly coupled through Compton interactions and 
collisions [50l. 
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fie flh 
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As for the fraction of free electrons, at the background level it obeys 
x' + 3HXr = Q , 



{AM) 



where the expression of Q can be found in Ref. [50| (more precisely what is called Q in 
this reference is equivalent to UcQ in here). It is in general a function of (xc, ric, T, H). 
T is the temperature of baryons and photons since it is approximately the same around 
recombination [50]. H is the local divergence of baryons {H = V^w{^/3) which at 
the background is the Hubble factor (in cosmic time). Q is obtained by taking the 
background value of all its arguments, that is 

Q = Q{xo,n,,f,H). (A.65) 

We report on Fig. lAll the numerical results obtained for the background fraction of free 
electrons. This has to be compared to the results obtained with RECFAST [96l |97] 
since we have implemented the correction to the three-level atom approximation 
implemented in this code, that is using its fudge factor of 1.14 for the last stages 
of the Hydrogen recombination and the fudge factor 0.86 for the Hel recombination. 
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Figure Al. Evolution of tlie free electrons fraction as a function of the redshift. 
We can see the transition of Hell and Hel at approximately z=6000 and z=2500 
respectively, and then the recombination of Hydrogen around z=1000. 



At first order, the evolution of the fraction of free electrons reads 



Xq 



simply 
(A.66) 



The perturbation Q^^'> is obained by taking the first order of Q{xc,nc,T, H) which 
is obtained by partial differentiation with respect to all the arguments of Q. The 
quantity H^^^^ is given by 




(A.67) 



The enhancement due to the perturbations in the fraction of free electrons is plotted 
in Fig. 
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Figure A2. Enhancement of the interaction rate due to the perturbations in the 
fraction of free electrons. S^i, Sy^ and 5r are respectively the relative perturbation 
of r', rie and pr- We plot 5^/ in continuous line, and S^, in dashed line. According 
to Eq. l lA.621 1. the difference between the two is the enhancement due to perturbed 
recombination. This has to be compared to the large scale limit of the delayed 
recombination S^,[l — l/3(lnXc)'/W], which is plotted in dashed-dotted line. We 
also plot in dotted line 3<5r/4 to show the departure from tight coupling around 
recombination. We have considered the two cases k = lO^^Mpc^^ (right) and 
k = 0.04Mpc-i (left). 



Appendix B. Emitting sources 



In the line of sight approach, we need to define the emitting sources. Instead of working 
with the brightness I and its corresponding electric and magnetic type polarizations 
£ and B, we work with 9 (defined in Eq. gU), £ = £/{4I) and B = B/{4I). They 
are expanded similarly to Eq. (|A.2ip . We obtain at first order 

1 



sg)(fc)^icW[£:]°(fc). 



(B.l) 



(B.2) 



However, if we want to have all effects but the late ISW localized on the last scattering 
surface, we need to integrate by parts in the line of sight integral solution. This boils 
down to using instead 

Sl^'^^ik) = b1 [*'(fc) + $'(fc)] + J?f'$(fc) [C(i)[l]?(fc) + T"'I°^'^(fc)] . (B.3) 
At second order the sources are 

SZfik) ^ <5^°, [$'(2) + *'(^)(fc)l + -5°<5«„f'$(^)(fc) (B.4) 
k_ 

1 



V3 

C(2)[Z]™(fc) 



_,^rn(2) 



(fc) 



s 



^(2) 



£,e 



S' 



m{2) 



ik) 



(k) 



C<-^^[£]T{k)+f'£f^\k) 
C^^'>[B]T{k) + f'B"'^^\k) 



(B.5) 
(B.6) 



The CMB bispectrum from the cosmological perturbations 41 

^ ic(i)(i)[i]r(fc) - , (B.7) 

^™(i)(i)(fc) = ic(i)(i)[£]r(fc) - , (B.8) 

s'i^'^^'Hk) ^ ic(i)(i)[6]r(fc) - i£(i)(i)[s]r(fc) . (B.9) 

Appendix C. Quadrupole initial conditions 

Since "H. ^ 1/77 in the radiation era, wc need to determine the initial conditions for 
up to terms in (A:?7i„it)'^, and thus to know the initial conditions for A/"™ up to terms 
of order kr\\^\\.- In fact, deep in the radiation dominated era, Nf'^^ ^ {kvY- 

= (c.i) 

ir(fc)=<^H^') = |$(^^(fc)-^{|$W(ki)$W(k,)} . (C.2) 

Using this, we deduce from the Boltzmann equation that at first order |70) 

<'\k)^^<^i^K (C.3) 
At second order we obtain the following initial condition 

^^''\k) = ^^('\k) + S, (C.4) 



(C.5) 



with 

- E 2i|^AAr<^^(k2)[5fc(")$(i)(ki) + (fcr^+fci"))*(i)(ki; 

n— — 1 J 

As for radiation, though the anisotropic stress is vanishing at initial time due to the 
high collision rate, the moment for £ = 2 is not vanishing since the bulk motion of 
radiation has a quadrupole [531 |67l |92] ■ We deduce that at initial time 

"'^jaS^ ^ S4Jr2i:i'tJ"«'"(ti)»'"('=2)| ■ (C.6) 
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Appendix D. Closure relations for vector and tensor modes 



We reproduce here the closure relations needed to integrate the Boltzmann equation for 
the vector and tensor modes. They are derived in Eqs. (5.4.48) and (5.4.49) of Ref. [80] 
which is not easily available. They correspond to the equations satisfied by the last 
multipoles in brightness and polarization kept in the truncated Boltzmann hierarchy. 
At first order we need only to retain the scalar closure relation {m = 0). However at 
second order we have to consider also the vector and tensor parts (m — 1, 2). Note that 
the presence of quadratic terms at second order implies that these closure relations 
arc not necessarily accurate as they were derived for the linear Boltzmann equation. 
For the brightness part, the closure relation reads 



ir' = k 



U+\m\ (2^+1), 
£-\m\ (21 -ly 
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iT" 
--Lc 



(D.l) 



For the polarization part, the closure relation for £™ and B™ is deduced by taking 
the real and imaginary part of 



(ff + \BT)' - k 
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m2 /^ + 2 (2^+1) 



J ( cm 

1) ^^^-1 



e-i) 



£ + 3 
krj 



{£]r + iSD 



(D.2) 



References 



[1 
[2: 
[3: 

[4 

[6; 

[7 

[s; 

[9 

[10: 
[11 
[12: 
[13: 

[14 

[is: 
[16: 

[17: 

[is: 

[19 

[20: 

[21 
[22: 

[23: 

[24: 



E. Kom atsu et al, (2010), |[arXlv: 1001 .4538]] 
Planck, 'http: //www.rssd. esa. int/ ' 

W. Hu and N. Sugiyama, Astrophys. J. 444, 489 (1995) , [JarXlv : astro-ph/9407093T| 

W. Hu and N. Sugiyama, Phys. Rev. D 51, 2599 (1995), [arXiv: astro-ph/9 411008T[ 

V. Mukhanov, H. Fcldman, and R. Brandcnbcrgcr, Phys. Rep. 215, 203 (1992^ 

A. Linde, Lect. Notes Phys. 738, 1 (2008), |[arXiv : 0705 . 0164] | 

S. Dodelson, Modern cosmology (2003). 

V. Mukhanov, Physical Foundations of Cosmology (2005). 

F. Bernardeau, Cosmologie (EDP Science, Paris, France, 2008). 

P. Peter and J. -P. Uzan, Primordial Cosmology (Oxford University Press, Oxford, UK, 2009). 
R. Durrer, The Cosmic Microwave Background (Cambridge University Press, 2008). 

S. Weinberg, Cosmology (Oxford University Press, 2008)^ 

J. Maldacena, JHEP 0305, 013 (2003), [ arXiv:astro-ph/0210603] [ 

M. AUshahiha, E. Silverstein, and D^Tong, Phys. Rev. D 70, 123505 (2004), 

[arXiv : hep-th/0404084] 
F. Bernardeau and J.-P. Uzan, Phys. Rev. D 66, 103506 (2002), [arXiv:hep-ph/0207295] 
E. Komatsu, The Pursuit of Non-Gauss ian Fluctuations in the Cos mic Microwave Background, 
PhD thesis, Tohoku University, 2002, | [arXlv : astro-ph/0206039] 

Bartolo, S. Matarres e, and A. Riotto, Phys. Rev. D 65, 103505 (2002), 
I [arXiv : hep-ph/01 1226lT| 

Bernardeau and J.-P. Uzan, Phys. Rev. D 67, 121301 (2003) , I [arXlv : astr o-ph/0209330T| 
Bernardeau, T. Brunier, and J.-P. Uzan, (2006), [arXiv : astro-ph/0604200^ 
I. Rigopoulos, E. P. S. She llard, and B. J. W. van Tent, Phys. Rev. D 73, 083522 (2006), 
[arXlv : astro-ph/0506704] | 

Sasaki, Progress of Theoretical Physics 120, 159 (2008) , ' [arXiv : 080 5 ■ 0974] | 
R. Bond, A. V. Prolov, Z. Huang, and L. Kofman, Physical Review Letters 103, 071301 



(2009) , [[arXi V : 0903 ■ 3407] 



N. Bartolo, E. Komatsu, S. Matarrese, and A. Riotto, Phys. Rep. 402, 103 (2004), 
I [arXi V : astr o-ph/0406398] 

E. Komatsu et al., Non-Gaussianity as a Probe of the Physics of the Primordial Universe and 
the Astrophysics of the Low Rcdshift Universe, in The Astronomy and Astrophysics Decadal 
Survey, Astronomy Vol. 2010, pp. 158— |-, 2009, [arXiv : 0902 . 4759] 



The CMB bispectrum from the cosmological perturbations 



43 



[25] N. Bartolo, S. Matarrese, and A. Riotto, (2010), |[arXiv: 1001 . 3957] | 

[26] N. Bartolo, S. Matarrese, and A. Riotto, JHEP 0404, 006 (2004), [arXiv : astro-ph/030 8088] | 
[27] T. Pyne and S. M. Carroll, Phys. Rev. D 53, 2920 (1996), [arXiv: astro-ph/9510041] 
[28] S. MoUerach and S. Matarrese, Phys. Rev. D 56, 4494 (1997), [arXiv : astro-ph/9702234] 
[29] D. M. Goldberg and D. N. Spergel, Phys. Rev. D 59, 103002 (1999), [arXiv:astro-ph/981125iri 
[30] E. Komatsu and D. N. Spergel, Phys. Rev. D 63, 063002 (2001), [arXiv : astro-ph/0005036] 
[31] D Babich, P. CremincUi, and M. Zaldarriaga, JCAP 0408, 009 (2004), 

[[arXlv : astro-ph/0 405356] 
[32] N. Bartolo, S. M atarrese, and A. Riotto, Phys. Rev. Lett. 93, 231301 (2004), 

[[arXlv : astro-ph/0 407505] 
[33] E. Komatsu, D. N. Sperge l, and B. D. Wandelt, Astrophys. J. 634, 14 (2005), 

[iarXlv:astro-ph/0305189] | 
[34] M. Liguor i, F. K. Hansen, E. Komatsu , S. Matarrese, and A. Riotto, Phys. Rev. D 73, 043505 
(2006), [iar"Xiv:astro-ph/0S09098T| ________ 

[35] N. Bartolo and A. Riotto, JCAP 0903, 017 (2009), [arXiv r08ri74584] 

[36] R. Khatri and B. D. Wandelt, Phys. Rev. D79, 023501 (2009), [arXiv : 0810. 4 370] | 

[37] L. Scnatore, S. Tassev, and M. Zaldarriaga, JCAP 0909, 038 (2009), [arXiv: 0812 .3658] ' 

[38] D. Nitta, E. Komatsu, N. Bartolo, S. Matarrese, and A. Riotto, JCAP 0905, 014 (2009), 

|[arXiv:0903. 0894] 

[39] D. Hanson, K. M. Smith, A. Challinor, and M. Liguori, Phys. Rev. D80, 083004 (2009), 
[arXiv: 0905. 4732] 

[40] L. Boubckcur, P. CremincUi, G. D'Amico, J. Norena, and F. Vernizzi, JCAP 0908, 029 (2009), 

[arXiv: 0906. 0980] 

[41] A. Mangilli and L. Verde, Phys. Rev. D80, 123007 (2009), [[arXiv: 0906 . 2317] I 

[42] M. Liguori, E. Sefusatti, J. R. Fergusson, and E. P. S. She llard, (2010), '[arXiv: 1001 . 4707] | 

[43] T. Okamoto and W. Hu, Phys. Rev. D66, 063008 (2002), [iarXlv: astro -ph/0206155] 

[44] N. Kogo and E. Komatsu, Phys. Rev. D73, 083007 (2006)7 fTar'Xiv : astr o-ph/0602099] 

[45] F. Bernardeau, T. Brunier , and J.-P. Uzan, Phys. Rev. D 69, 063520 (2004), 

|[arXlv:astro-ph/0311422] | 
[46] S. Weinberg, Phys. Rev. D 72, 043514 (2005), [arXiv:hep-th/0506236] 
[47] S. Weinberg, Phys. Rev. D 74, 023508 (2006), [arXiv:hep-th/0605244]' 
[48] S. Weinberg, Phys. Rev. D 78, 063534 (2008), [arXiv : 0805 . 3781] 

[49] C. Pitrou, J.-P. Uzan, and F. Bernardeau, Phys. Rev. D 78, 063526 (2 008),|[arXlv:0807.0341] [ 
[50] L. Senatore, S. Tassev, and M. Zaldarriaga, JCAP 0908, 031 (2009), [arXiv: 0812 . 3652]J 
[51] R. Khatri and B. D. Wandelt, Phys. Rev. D79, 023501 (2009), [arXiv: 0810 . 4 370] 
[52] K. Nakamura, Prog. Theor. Phys. 117, 17 (2007), [arXiv : gr-qc/0605108] 
[53] C. Pitrou, Class. Quant. Grav. 26, 065006 (2009), [arXiv: 0809 . 3036] 
[54] C. Pitrou, Gen. Rel. Grav. 41, 2587 (2009), [arXiv: 0809 . 3245] 

[55] N. Bartolo, S. Matarrese, and A. Riotto, JCAP 0508, 010 (2005), [arXlv:astro-ph/0506410] [ 
[56] N. Aghanim, S. Majumdar, and J. Silk, Reports on Progress in Physics 71, 066902 (2008), 
JarXiv:07 11.0518] 

[57] C. Pitrou, F. Bernar deau, and J.-P. Uzan, (2009), [arXiv : arX iv : 0912 ■ 3655] | 
[58] C. Pitrou, [http: //Tcg.port . ac .uk/~pitrouc/cmbquick.htm 

[59] M Bruni, S. M atarrese, S. Mollerach, and S. Sonego^ Class. Quant. Grav. 14, 2585 (1997), 

I [arXiv : gr-qc/9609040] 
[60] WMAP, E. Komatsu et al, Astrophys. J. Suppl. 180, 330 (2009), [arXiv : 0803 . 05 47]] 
[61] J. Lesgourgues and S. Pastor, Phys. Rcpt. 429, 307 (2006), [arXiv:astr o-ph/0603494] | 

[62] I. Tereno et al, Astron. Astrophys. 500, 657 (2009), [arXiv: 0810 . 0555] | 

[63] W. Hu and M. White, Phys. Rev. D 56, 596 (1997), ,[arXiv: astro-ph/9702170] | 

[64] A Coc, K. A. Olive , J.-P. Uzan, and E. Vangioni, Phys. Rev. D79, 103512 (2009), 

|[arXlv:0 811.1845] 

[65] M. Bucher, K. Moodley, and N. Turok, (2000) , [[arXlv : astro-ph/001 1025T| 

[66] A. Linde and V. Mukhanov, Phys. Rev. D 56, 535 (1997), [arXlv: astro-ph/9610219] | 
[67] C. Pitrou, Class. Quant. Grav. 24, 6127 (2007), [arXiv: 0706 . 4383] 1 

[68] K. A. Malik and D. Wands, Class. Quant. Grav. 21, L65 (2004), [arXiv: astro-ph/0307085] | 

[69] F. Vernizzi, Phys. Rev. D 71, 061301 (2005), [arXiv:astro-ph/0411463] I 

[70] C.-P. Ma and E. Bertschinger, Astrophys. J. 455, 7 (1995), [arXiv : astro-ph/9506072] 

[71] P. CremincUi, A. Nicolis, L. Senatore, M. Tcgmark, and M. Zaldarriaga, JCAP 0605, 004 

(2006), [arXlv:astro-ph/0509029] 
[72] P. CremincUi, L. Scnatore,_M. Zaldarriaga, and M. Tegmark, JCAP 0703, 005 (2007), 
[[arXiv : astro -ph/0610600] | 



The CMB bispectrum from the cosmological perturbations 



44 



[73] M. Zaldarriaga and U. Scljak, Phys. Rev. D 55, 1830 (1997), |[arXlv: astro- ph/9609170T| 

[74] W. Hu and M. White, New Astron. 2, 323 (1997), [arXivrastro-ph/9706147r] 

[75] W . Hu, U. Se ljak, M. J. White, and M. Zaldarriaga, Phys. Rev. D 57, 3290 (1998), 

|[arXlv:astr'o -ph/9709066] 
[76] J. R. Bond, Theory and Observations of the Cosmic Background Radiation, in Cosmology and 

Large Scale Structure, edited by R. Schaeffer, J. Silk, M. Spiro, & J. Zinn-Justin, pp. 469, 

1996. 

[77] C. Pitrou, J. -P. Uzan, and F. Bernardea u, in p reparation. 
[78] L. R. Abramo and T. S. Pereira, (2010), [iarxTv : 1002 .3173] | 

[79] S. Prunet, J.-P. U zan, F. Be rnardeau, and T. Brunier, Phys. Rev. D71, 083508 (2005), 

JarXiv:astro -ph/0406364] I 
[80] A. Riazuelo, Signature de divers modeles d'univers primordial dans les anisotropics du 

rayonncment fossile, PhD thesis, Universitc Paris Sud, 2000. 
[81] E. Fenu, C. Pitrou, and R. Maartens, in preparation. 
[82] R. Khatri and B. D. Wandelt, (2009), [arXiv : 0903 . 0871] | 

[83] M Beneke and C. Fidler, (2010), [arXiv : 1003 . 1834] |. 

[84] W. Hu, Phys. Rev. D 62, 043007 (2000), [arXiv: astro-ph/0001303] ' 

[85] N. Bartolo, S. Matarrese, and A. Riotto, JCAP 0606, 024 (2006), [arXiv:as tro-ph/0604416] [ 

[86] M. Zaldarriaga, Phys. Rev. D62, 063510 (2000), [arXiv : astro-ph/9910498] ' 

[87] A. Lewis and A. Challinor, Phys. Rept. 429, 1 (2006), [arXiv : astro-ph/0601594] I 

[88] D. N. Spergel and D. M. Goldberg, Phys. Rev. D 59, 103001 (1999), [arXiv:astro-ph/981 1252] | 
[89] L. Verde and D. N. Spergel, Phys. Rev. D 65, 043007 (2002), [arXiv : astro-ph/0108179f] 
[90] W. Hu and A. Cooray, Phys. Rev. D63, 023504 (2001), |'[arXlv : astro-ph/0008001] j 
[91] V. Acquaviva, N. Bartolo, S . Matarrese, and A. Riotto, Nucl. Phys. B 667, 119 (2003), 

|[arXiv:astro-ph/0209156] | 
[92] N. Bartolo, S. Matarrese, and A. Riotto, JCAP 0701, 019 (2007) , j [arXlv : astro-ph/06101 loT| 

[93] P. J. E. Peebles and J. T. Yu, Astrophys. J. 162, 815 (1970)^ 

[94] A. Lewis and A. Challinor, Phys. Rev. D76, 083005 (2007) , ' [arXiv : astro-ph/0702600] ' 
[95] B. Novosyadlyj, Mon. Not. Roy. Astron. Soc. 370, 1771 (2006), [arXiv:astro-ph/0603674] 
[96] S. Seager , D. D. Sassel ov, and D. Scott, Astrophys. J. 523, LI (1999), 

[arXiv:a s"tro-ph/9909275] | 
[97] W. Y. Wong, A. Moss, and D. Scott, (2007), |[arXiv: 0711 ■ 1357] | 



